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We systematically study the unification of gauge couplings in the presence of (one or more) effec- 
tive dimension-5 operators cHG^vG^" /AMpi, induced into the grand unified theory by gravitational 
interactions at the Planck scale. These operators alter the usual condition for gauge coupling uni- 
fication, which can, depending on the Higgs content H and vacuum expectation value, result in 
unification at scales Mx significantly different than naively expected. We find non-supersymmetric 
models of SU(5) and 5*0(10) unification, with natural Wilson coefficients c, that easily satisfy the 
constraints from proton decay. Furthermore, gauge coupling unification at scales as high as the 
Planck scale seems feasible, possibly hinting at simultaneous unification of gauge and gravitational 
interactions. In the Appendix we work out the group theoretical aspects of this scenario for SU(5) 
and 50(10) unified groups in detail; this material is also relevant in the analysis of non-universal 
gaugino masses obtained from supergravity. 
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I. INTRODUCTION AND OUTLINE 

What are the boundary conditions for grand unification? One typically assumes that the gauge couplings of the 
broken subgroups must become numerically equal at the unification scale Mx PQ. However, effects from physics 
above the unification scale can alter the gauge coupling unification condition. In an effective field theory approach, 
such effects can be caused by dimension-5 operators of the form cHG / AMpi, which shift the coefficients of the 
gauge kinetic terms in the low-energy theory after the Higgs H acquires a vacuum expectation value in grand unified 
symmetry breaking [21 [3]; one obvious source of such operators is quantum gravitational effects. Indeed, it would be 
unnatural (or require some special explanation) to assume that the Wilson coefficients c above be zero or especially 
small 0]; the default assumption should be that these coefficients are of order unity in grand unified models, with 
consequent unification conditions. 

In conventional unification models, one might expect (H) ~ 10 16 GeV, plausibly leading to effects from quantum 
gravity of order a fraction of a percent, (H)/Mpi ~ 10~ 3 , on the gauge coupling unification condition. In [5] we 
showed that these dimension-5 operators can be even more relevant than previously suspected since the Planck mass 
Mpi tends to be smaller than naively assumed due to its renormalization group evolution (6l [7] under the influence of 
the large number of fields in supersymmetric grand unified theories. It was noted [5] that these dimension-5 operators 
introduce in supersymmetric unification models an uncertainty that can be bigger than the two-loop effects which are 
considered to be necessary to obtain good numerical unification of the gauge couplings. 

The aim of this paper is different. We study whether the dimension-5 operators discussed above can lead to perfect 
gauge coupling unification without supersymmetry by their modifying of the gauge coupling unification condition. 
This unification scheme has been studied previously in the literature for models with and without supersymmetry, 
e.g. in [2l 151 [51 15HT5] . but in less detail and generality and mostly only the effect from a single gravitational operator 
has been considered. 

In particular, in this paper we examine, in a systematic way, the effects when two or more dimension-5 opera- 
tors are present in a theory. Unification under multiple dimension-5 operators has been studied before [13] for the 
supersymmetric case |16j . for which, however, viable gauge coupling unification is well known |17H19] . The main 
result of the present paper is that the measurement of the gauge couplings at the Z mass [20] is compatible with 
non-supersymmetric grand unification based on 517(5) or bigger groups like 5(7(10). That is, grand unification does 
not require a supersymmetric extension of the standard model for a range of natural values of the Wilson coefficients 
c. 
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As opposed to models with one dimension-5 operator, in unified models with two or more such operators, the 
unification scale Afj can be varied in a continuous and controlled manner with the Wilson coefficients c and Higgs 
vacuum expectation values (vevs), and we examine this quantitatively in SU(S) and S'O(IO) theories, primarily for 
the non-supersymmetric case (cf. [13) for the supersymmetric case). We find that unification under this scheme can 
naturally raise the unification scale Mx much above the conventional 10 15 GeV to 10 16 GeV, and even allows to obtain 
gauge coupling unification at or near the Planck mass which is suggestive of unification of all gauge and gravitational 
forces at a common scale (see also [2T][22]). Any such unification is safe from the proton decay constraint. This, and 
the fact that exact gauge coupling unification in the first place can be accomplished, may be a useful tool for model 
building. 

Note, due to the gravitational origin of these dimension-5 operators, the unification scheme considered here is distinct 
from the Lavoura-Wolfenstein result [33J or similar effects according to which non-supersymmetric unification 
can result from particle thresholds. In particular, as can be seen from the form of these effective operators, their 
effect increases as the unification scale gets closer to the Planck scale, i.e., as the Higgs vevs get bigger, which in turn 
allows to self-consistently shift Mx to values much larger than conventionally assumed. 

Realistic supersymmetric unification theories have issues linked to their large particle content. One is the presence of 
Landau poles between the unification and the Planck scale, another is possible lack of calculability in these theories due 
to potentially large running of the Planck mass [5]. Both issues are avoided in their non-supersymmetric counterparts. 
On the other hand, if supersymmetry is abandoned, the unification scale needs to be stabilized with respect to the 
gravitational scale (as, of course, does the weak scale). This "little hierarchy problem" can be avoided by one of the 
scenarios mentioned above where grand unification happens close to the Planck mass, thereby reducing or eliminating 
the hierarchy between the two scales. 



Outline of the paper. For a physical description of the scenario and effects and for exemplary numerical results, see 
the Introduction and Summary, Sections [I] and VIII More detail is given in the main text, and the group theoretical 
formalism in the Appendix. 

In SectionlllJ we introduce the dimension-5 operators cHG /4Mpi under investigation, first specifically for the 
case of an SU[5) grand unified theory, and also discuss possible sensible choices for the Planck scale Mpi suppressing 
these operators. We then describe how these operators modify the usual gauge coupling unification condition. This 
section sets up the notation necessary to understand the figures, tables and most of the details in the main text. 

In Section pTI} we look at numerical results for possible gauge coupling unification under the modified gauge coupling 
unification condition in non-supersymmetric SU(5) models. First (Section III A) we review the known effects in the 
presence of a single such operator [5J [3J [U 115) , also addressing uncertainties in low-energy measurements of the gauge 
couplings and modifications due to two-loop running. Then, in Section III B it is numerically shown that (and how) 



the unification scale may be continuously varied in models with two different dimcnsion-5 operators. This, being one 
of the main results of this paper, is also true in models with more such operators, and we give the general treatment 
in Section III C also deriving an estimate to assess which unification scales Mx are achievable naturally. 

With these numerical results and estimates, in Section [TV| we see that even the lowest unification scales Mx, that 
are achievable in a natural way through our effect, are not in conflict with the current lower bounds on the proton 
lifetime, coming from the non-observation of proton decay so far. Thus, the dimension-5 operators can facilitate exact 
gauge coupling unification without supersymmetry and also allow the proton decay limit to be evaded (unlike in naive 
SU(5) unification). We also investigate how further improvements in the proton lifetime bound constrain our models. 

In Section |Vj we note that numerical gauge coupling unification is possible near or at the Planck scale in a very 
natural way in models with two or more dimcnsion-5 operators. We speculate that this might hint at unification of 
the gauge interactions and gravity at a common scale ("gauge-gravity unification", see [21]). 

Section VI contains the scenario for an S'O(IO) grand unified gauge group. The setup and unification results for 
exemplary cases are described, analogous to the SU (5) case in the preceding sections. The main differences to SU(5) 
are emphasized, namely the fact that a continuously varying Mx can be achieved with merely a single dimension-5 
operator, and the possible role of these operators in 50(10) breaking with intermediate scales is briefly described. 

In Section |VII| we briefly look at the effect of the dimension-5 operators in models of supersymmetric unification 
[TBHTS] . making contact to the literature (e.g., [5] [T3J [T5] ) • We find that in supersymmetric models the unification 
scale can also be shifted around, although unification close to the expected 2 x 10 16 GeV seems most likely. The main 
part of the present paper focuses on non-supersymmetric models as it is those where viable unification is commonly 
believed to be difficult or impossible (for further reasons, see the last paragraph of Section pi]) . 

We conclude and summarize the main ideas and mechanisms in Section |VIII[ illustrated with some numerical 
examples showing the size of the effects. 

In an extensive Appendix, we present all normalization conventions and group theoretical details. In particular we 
carefully derive and give all of the Clebsch-Gordan coefficients $? r ) t associated with SU(5) and 50(10) breaking (the 
latter case takes up most of the work) to the standard model, in different bases {t}; only the coefficients relating to 
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the standard model gauge group factors s have been given before in the literature, see especially [53], and often only 
their relative sizes, whereas we here follow a uniform absolute normalization scheme across different representations 
r. These analytical results are also useful for the treatment of non-universal gaugino masses obtained from M = 1 
supergravity models of 577(5) or 5*0(10), for which the group theory involved is very similar, see, e.g., (551 123- The 
Appendix furthermore describes the modified gauge coupling unification condition and a systematic method (which 
was used for the numerical work in this paper) for its solution. 



II. DIMENSION-5 OPERATORS AND UNIFICATION CONDITION - SU(5) CASE 

In this section we describe, specifically for the case of non-supersymmetric SU (5) grand unification models [28] , 
first the operators necessary for our scenario and then their effect on the condition for gauge coupling unification. The 



50(10) case holds additional subtleties and possibilities, to be deferred to Section VI We present actual numerical 



unification results for the SU(5) case in Section III and elaborate on physically interesting scenarios in Sections IV 
and El 

At the energies available to us in present-day particle physics experiments, Nature is described very well by non- 
Abelian gauge field theories which are non-gravitational, although we know that at some energy scale a more complete 
description must take gravity into account. In other words, at our energies the effective Lagrangian of the world is 
a gauge theory with certain additional non-renormalizable operators of dimension 5 and higher, accounting for the 
gravitational interactions that have been integrated out. These operators must obey the symmetries (gauge and 
Lorentz invariance, etc.) of the low-energy theory and are suppressed by powers of the gravitational scale Mpi, so 
they are seemingly negligible at our energies. By this logic, grand unification also appears as an effective theory, valid 
between the unification scale Mx and well below the Planck scale Mpi, and contains higher-dimensional operators 
induced by gravity and suppressed by Mpi; at energies ~ Mx, however, such operators are potentially much more 
significant due to the proximity of the scales Mx and Mp; . 

One set of such dimension-5 operators, that may have important effects in grand unification, are singlets formed 
from gauge field strengths G^ and Higgs multiplets Hi of the grand unified gauge group G = SU(5), 

C = -r^H? b G*G b ^ , (1) 

suppressed by one power of the Planck mass Mpi such that Cj are dimensionless (Wilson) coefficients. The index i 
is summed implicitly and includes the possibility that the effective Lagrangian may contain several such operators 
involving different Higgs multiplets Hi of the theory, which will be one of our main tools later on. In an SU (5) gauge 
theory, operators can be formed gauge-invariantly only with Hi in the representations r; = 1, 24, 75 or 200 
(although the theory might contain additional multiplets in other representations); these irreducible representations 
(irreps) can all uniformly be written in component notation with symmetric adjoint indices a, b of the gauge group, 
which establishes a common normalization for the different operators i in ([I]). (Here and later, see the Appendix for 
a careful treatment of the relevant group theoretical aspects and normalization conventions; see also |29|.) 

The dimension-5 operators ([I]) are suppressed by the Planck scale Mp;, the energy scale at which quantum gravity 
sets in, which, as an interaction not accounted for by the renormalizable terms in the Lagrangian, induces these 
effective operators. There is some ambiguity (or arbitrariness) as to whether the appropriate suppression scale Mpi 
should be taken to be the "naive" Planck scale G JV 1 '' 2 = 1.2 x 10 19 GeV or, more commonly, the reduced Planck scale 

(87rGAr) -1 ^ 2 = 2.4 x 10 18 GeV, as this is the quantity that controls quantum gravity computations. To leave this 
choice open, we parametrize 

Mm m _ 1.2,10-G.V . (2) 

£ = 1 corresponds to the choice of the naive Planck scale as the suppression scale, £ rc d = \/8tt m 5 to the reduced 
Planck scale, which is what we assume (implicitly) in most discussions. The smaller a suppression scale Mp\ one 
accepts, i.e. the bigger £ one chooses, the more pronounced the effects of the operators will be, at fixed Wilson 
coefficients Ci . Effects of equal size can be achieved for indirectly proportional coefficients a — > Ci /£ when changing 
£. (Conservatively, we put an explicit factor 1/4 in ([I]) to avoid overcounting of terms in the contraction of two gauge 
field strengths, just as in the canonical gauge boson kinetic term.) 

Also concerning the choice of an appropriate suppression scale Mpi, we have shown previously [51 [7] that the 
fundamental value of Newton's constant (i.e., at high energies) is different from its observed low-energy value Gn 
used in d2h : matter field fluctuations of N real scalar, Ni/ 2 Weyl fermion and JVi gauge boson fields lead to a running 
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of Newton's constant 



GQj.) 



1 

Gn 



12tt 



(iVo + N 1/2 - 4/v\) 



(3) 



at one loop, similar to the running of gauge couplings (see also [51 1301 |3"T] ) . Then determining the fundamental gravita- 
tional scale Mpi via G(/z = Mpi)^ 1 ^ 2 ~ Mpi yields a value lower by a factor of 77 = ^/l + (Nq + Ni/ 2 — 4/Vi)/12-7r than 
without the running. This additional change in the suppression scale Mpi is easily incorporated into our parametriza- 
tion ([2| via £ — > £ run = £77, and Table [I] illustrates the size of this running effect in various grand unified models. 
Furthermore, as illustrated in the last row of the table, appropriate choices of £ can accommodate suppression scales 
Mpi of other origin as well, e.g. string compactification scales as in [3J. 



particle physics model 


jV = N + N 1/2 - 47Vi 


r, = y/1 + 7V/127T = f un 


= V^V 


log 10 M„p/Mx 


no running of Gat 




1 


5.0 




standard model 


1 


1.0 


5.1 




5(7(5) w/ 5, 24 


-17 


0.74 


3.71 




577(5) w/ 5, 200 


159 


2.3 


11.5 


42 


577(5) w/ 5, 24, 75 


58 


1.6 


8.0 




577(5) w/ 5, 24, 75, 200 


258 


2.8 


14.0 


14 


5O(10) w/ 10, 16, 45 


-35 


0.27 


1.34 




5O(10) w/ 10, 16, 210 


130 


2.1 


10.6 




5O(10) w/ 10, 16, 770 


690 


4.4 


22.0 


3.9 


SUSY-5C/(5) w/ 5, 5, 24 


165 


2.3 


11.6 




SUSY-5(7(5) w/ 5, 5, 24, 75 


390 


3.4 


16.9 


3.6 


SUSY-5!7(5) w/ 5, 5, 200 


693 


4.4 


22.1 


0.85 


SUSY-SO(IO) w/ 10, 16, 16, 45, 54 


432 


3.5 


17.7 


11 


SUSY-5O(10) w/ 10, 16, 16, 210 


765 


4.6 


23.1 


1.8 


SUSY-5O(10) w/ 10, 16, 16, 770 


2445 


8.1 


40.7 


0.27 


compactification scale M c (e.g. [3]) 




up to £ ~ 100 





TABLE I: For models with various particle contents, the third column shows the effect Mpi — > Mpi/rj on the Planck scale 
entailed by the running |3| of Newton's constant. The fourth column gives numerical values for £ in |2| if the additional 
reduction factor V87F from perturbative quantum gravity is taken into the suppression scale, as commonly done. Throughout 
we assume three generations of fermions, and grand unified models are characterized by their gauge group and Higgs content. 
In almost all cases - and certainly in all models of our main interest, namely in non-supersymmetric grand unified models 
with several different Higgs multiplets at the unification scale - reasonable values for Mpi are smaller than the naive value 
1.2 xlO 19 GeV, in some cases by as much as roughly an order of magnitude, i.e. 1 < £ < O(10) in The last column quantifies 
roughly how many orders of magnitude above the unification scale Mx the grand unified gauge theories enter non-perturbative 
strong coupling regime gc(M np ) = i/4-k, a dash ( — ) indicating that the theory is asymptotically free; here, we have assumed 
go(Mx) 2 /4-7T = 1/40 for the non-supersymmetric and (conservatively) ga(Mx) 2 /4-7T = 1/30 for the supersymmetric cases. 
When M np > Mpi, the unified field theory has a chance of describing Nature perturbatively up to the onset of quantum gravity 
at Mpi. 

Because they respect the symmetries of the theory and because quantum gravity effects are mediated at the scale 
Mpi, operators ([lj should be expected in the effective Lagrangian of grand unification and, although their sizes cannot 
be computed in this effective field theory point of view, they are expected a priori with natural Wilson coefficients of 
order |cj| ~ 0(1) (for values much bigger than this, the proposed effective theory is not a good low-energy description 
and one loses perturbative control, whereas |cj| <C 1 seemingly constitutes fine-tuning [3]). And even though irrelevant 
for our purposes here, there are several known ways to generate these operators. For example, they arise from an J\f = 1 
supergravity ultraviolet completion with non-canonical gauge kinetic function f ab (Hi) as the lowest order (non-trivial) 
terms in an expansion of f ab in Mpi |26U27j (in this scenario, supersymmetry needs to be broken in the grand unified 
theory at Mpi to conform to our non-supersymmetric analysis). Gravitational instantons can induce such effects 
also |32j . And in descending from higher-dimensional completions to four dimensions, spontaneous compactification 
generates operators ([T]) suppressed by the compactification scale M c = Mpi [33] . At any rate, lacking knowledge of 
quantum gravity, it seems most reasonable to assume the presence of effective operators |l]) with coefficients Cj of 
order one. 
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We now describe the effect of the operators ([I]) on the condition for gauge coupling unification. At the scale Mx of 
grand unified symmetry breaking, some Higgs multiplets Hi acquire nonzero vacuum expectation values (vevs). For 
simplicity and definiteness we assume that all Higgs fields, except for the multiplet containing the standard model 
Higgs, acquire vevs at the scale Mx; as will become clear, other multiplets, that get nonzero vevs only at lower 
scales, contribute proportionally less, but can in principle be treated equivalently (this assumption also avoids the 
introduction of further mass hierarchies into the model). 

The vevs (Hi), acquired well above the electroweak scale, have to be invariant under the standard model subgroup 
G321 = SU(3)c x SU(2)l x U(1)y C SU(5). For each multiplet Hi that can occur in the dimension-5 operators 
|l]), this requirement determines its vevs (Hf h ) up to an overall scale Vi (see the Appendix). Replacing Hi by their 
vevs, the operators ([!]) modify the kinetic terms of the gauge bosons in the Lagrangian at the unification scale Mx 
by adding to them: 

C = -\Gl v G a ^ +Y (Hf) G a uv G b ^ 
4 Z— < AMpi M 

i 

= ~ \0- + - + ^Wuh + «i) VSp) + • • • , (4) 

where the ellipses denote the non-standard model gauge bosons of 577(5), which become massive and (quasi) non- 
dynamical below Mx- The e s indicate the modifications to the gauge kinetic terms, 

e '=Ej| L 7^ <) (for s = 3, 2,1); (5) 

i 

they generically differ for each factor s = 3, 2, 1 of the standard model gauge group G321 and depend on the Higgs 
content Hi and on the sizes CiVi/Mpi of the vevs and Wilson coefficients relative to the suppression scale (thus, a Higgs 
that acquires its vev v ~ Mj <C Mx at an intermediate scale contributes much less than a Higgs at Mx)- The Si are 
group theoretical factors (Clebsch-Gordan coefficients) specific to the embedding G321 C SU(5) and characterize the 
possible standard model singlet vevs (Hi); they depend only on the representation rj of Hi and are given in Table III] 



SU (5) irrep r 


5[ r) 


r(r) 
°2 


8i r) 


1 


-1/^/24 


-1/V24 


-V\ 


/24 


24 


1/V63 


3/V63 




/63 


75 


5/V72 


-3/V72 




/72 


200 


-10/\/l68 


-2/ v / 168 


-i/v 


168 



TABLE II: The Clebsch-Gordan coefficients 5 S in |5]l associated with the embedding G321 C SU(5), for each irrep r of 57/(5) 
that can occur as a multiplet Hi in (m). Sensible normalization conventions (see the Appendix), which ensure uniform treatment 
of different dimension-5 operators (fil, fix the values of all 5^ up to an overall sign for each r. 

In the effective Lagrangian below the unification scale Mx, one would like to have canonically normalized gauge 
fields as opposed to Q, since it is the coupling constants associated with those that obey the familiar renormalization 
group (RG) equations (/3-functions) and that are measured in low-energy experiments. This can be achieved by a 
finite (and usually small, see later) redefinition A?, — s- (l + e s ) 1//2 A^ of the gauge fields associated with each standard 

model factor s = 3, 2, 1, which has to be accompanied by a redefinition g s — > (1 + e s ) _1 / 2 g s , so as not to affect the 
interaction strength. Gauge coupling unification requires that at the unification scale Mx the couplings before this 
redefinition meet at a common value go = gc(Mx), the gauge coupling of the unified group G at the unification 
scale. In terms of the rescaled couplings a s = g 2 /47r, this reads: 

(l + ei)ai(Mjr) = (l + e 2 )a 2 (M x ) = (l + e 3 )a 3 (M x ) = ^ = a G . (6) 

Under our assumption of one-step breaking at Mx to the standard model, the running coupling functions a s (fi < Mx) 
are fixed by the fairly precise low-energy measure ments [2T)| and their RG evolution [T7]; the low-energy values and 
/3-functions to one loop are given in the Appendix ( A.30|-( A. 32 1. 



In Section III we will be looking for non-supersymmetric models of SU(5) grand unification (specified by their 



Higgs content Hi) which feature exact unification of the standard model gauge couplings a s (fi) under the unification 
condition (ml. Such gauge coupling unification can happen naturally at scales Mx > 10 16 GeV, larger than normally 
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expected, so as to escape the proton decay limit usually encountered in non-supersymmetric grand unification, and 
even at scales as large as the Planck scale. Through several criteria we will assess reasonableness of these models, in 
particular by checking whether their required Wilson coefficients have natural sizes |c,| ~ 0(1) and also by looking at 
the modifications to the gauge kinetic terms in @ (e.g., if e, < -1 for any s, the gauge kinetic terms in the broken 
theory would have the wrong sign). 

One other criterion for all models, necessary to sensibly claim unification, is the requirement that the masses of 
all non-standard model gauge bosons ("superheavy gauge bosons") be close to the scale of grand unified symmetry 
breaking, i.e., close to the unification scale Mx- This is to ensure rightful use of the standard model RG equations 
with three differently evolving gauge couplings up to the scale Mx , where the gauge couplings are to unify at go ((6} , 
and of the running of a single coupling of the unified gauge theory thenceforth with all gauge bosons being massless. 
At least this is necessary in the absence of intermediate scales, which is usually understood for SU (5) unification 
models (see above). In SU(5), all superheavy gauge bosons get equal masses from grand unified symmetry breaking, 




where the sum now runs over all Higgs multiplets i in the theory that acquire nonzero vev Vi at Mx, with the Casimir 
invariants Ci of the representations (cf. Table |V|) . For definiteness, and to conform to standard treatment, we take 
the above requirement then to be 

M gb = M x ■ (8) 

Furthermore, since Higgses in any of the representations 24, 75, 200 can achieve grand unified symmetry breaking 
SU(b) — > G321, we will make the simplifying assumption that the theory only contain Higgses able to form the 
dimension-5 operators ([I]); the more general case is not harder to treat but increases the particle content of the 
models and requires larger Wilson coefficients to achieve the desired effects. 

A requirement similar to Q should be put on the masses of the superheavy scalars as well, i.e. the Higgs fields 
that acquire nonzero vev at Mx, which we, however, will neglect since generically most superheavy scalars are more 
massive than the (lightest of the) superheavy gauge bosons and since the superheavy scalars do not influence the 
running of the standard model gauge couplings at one loop. Also, we neglect heavy particle thresholds [231 [H], to 
demonstrate unification due specifically to the dimension-5 operators ([!]) as a proof of principle, cleanly separated 
from other effects. 



In this paper we focus on non-supersymmetric grand unification. Contrary to the standard lore, and this is one of 
our main findings (Sections III and IV I, the altered unification condition ^ can yield successful non-supersymmetric 
unification models which satisfy the proton lifetime constraint. Furthermore, non-supersymmetric grand unified mod- 
els hold the attractive hope of describing physics up to the Planck scale in a perturbative way, whereas supersymmetric 
models commonly become strongly coupled before reaching the Planck scale. For several non-supersymmetric as well 
as supersymmetric grand unified models, the right column in Table [I] indicates roughly how many orders of magnitude 
above the unification scale Mx the theories become non-perturbative (i.e., near a Landau pole) and illustrates that 
generally only non-supersymmetric models are safe in that regard, requiring roughly 3 orders of magnitude between 
Mx and the non-perturbative regime M np . Furthermore, as they contain fewer scalars and fermions, the running (13J) 
of Newton's constant and the associated change of the Planck scale Mpi from its naive or reduced value tend to be 
smaller in the non-supersymmetric models (see Table |l| , diminishing uncertainties in these models [5] . 



III. MODELS AND UNIFICATION RESULTS IN THE NON-SUPERSYMMETRIC SU(5) CASE 

In this section, we will quantitatively examine the effects of multiple gravitationally induced dimension-5 operators 
([!]) and of the subsequently modified condition ^ for gauge coupling unification in several non-supersymmetric SU(5) 
models. All models are presumed to have the minimal fermionic content of three standard model families, but differ in 
their Higgs content responsible for grand unified symmetry breaking; without loss (see above) we only consider Higgs 
representations 1, 24, 75, 200. As the Higgs multiplets under consideration are suitable for breaking SU(5) down to 
the standard model, we do not consider the Higgs (scalar) potential in these unified theories explicitly; rather, we take 
as the parameters of the models the Higgs vevs Vi directly, which are acquired in grand unified symmetry breaking as 
a consequence of the Higgs potential, and the Wilson coefficients Cj of the dimension-5 operators 0. 

Our method is then as follows. After specifying a grand unified model by its Higgs content, we scan its parameter 
space {ci,Vi} for points that, below the breaking scale Mx, yield the actual running gauge couplings a s {n) of the 



7 



standard model. Phrased in a bottom-up language, we are looking for points {ci,Vi} in parameter space that result 
in a unification condition ([6]), according to which the actually observed gauge couplings a s (/j,) of the standard model 
unify, while simultaneously requirement ^ holds, namely that the superheavy gauge boson masses ^ be equal to 
the unification scale. In fact, any such point {cj,i>i} determines Mx and the unified gauge coupling go = gc(Mx) 
uniquely; sec the Appendix below (A. 33) where we also outline how all points {ci,Vi} may be found. 

In the next subsection we look at the known scenario in which the theory contains only one Higgs that can form 
a dimension-5 operator, as a warm-up and comparison to established results [21 131 [TB]. For the case of two or more 
operators 0, we find points in parameter space that yield physically viable unification. These effects are described in 
detail exemplarily for scenarios with two dimension-5 operators in Section III B before moving on to a more general 
treatment of models with any number of dimension-5 operators where we also pay special attention to naturalness. 



A. Unification in models with one dimension-5 operator 



When the modifications e s to the gauge kinetic terms Q come from only one Higgs Q, their ratio €\ : 62 ■ £3 is 
completely determined, see ([5]). For the given running gauge coupling functions a s (ii) of the standard model, this 
ratio uniquely determines, via the unification condition ([6]), the possible unification scale Mx, the unified coupling 
a G — OgI^ an d the required absolute sizes of all e s ; subsequently, the required Higgs vev v and Wilson coefficient 
c necessary for unification may be computed. It then has to be checked whether all these values are physically 
reasonable. 



For numerical results, see Table III which shows these quantities for the SU(5) models with a sole Higgs multiplet 
24, 75, or 200. A sole singlet 1 cannot modify condition Q to make the standard model gauge couplings unify 
since t\ — ti = e 3 in this case, neither can a singlet vev break the unified symmetry in the first place. Although 
the required modifications e s to the gauge kinetic terms are not problematic (in particular, e s > —1) and even small, 
the necessary Wilson coefficients c are big in the case of a 75 and especially large for a 24 Higgs, regardless of any 
reasonable choice of suppression scale; operators with such large coefficients are not expected naturally in an 
effective theory and might preclude a perturbativc treatment. Furthermore, in the case of a 24 the unification scale 
Mx would be far too low to satisfy current constraints on the proton lifetime, and possibly slightly too low for a 75 
(for details, see Section IV). For a 200 Higgs the Wilson coefficient c has natural size and, depending on convention 
([2]), the unification scale may still lie below or at the Planck scale Mpi, which may be taken as a hint to a scenario 
of simultaneous gauge-gravity unification (see Section [V| . 
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53.4 


0.53/5 


1.1 x 10 19 


0.363 


9.8 x 10 17 


2.6/5 


1 


5 


2.3 


11.5 


0.23 



TABLE III: For all suitable Higgs irreps, this table shows the parameters a non-supersymmetric 5(7(5) model must have if 
unification of the standard model happens by means of mechanism |6]) with only one dimension-5 operator ([I]). Columns 2-6 are 
obtained by using the one-loop /3-functions for the standard model gauge couplings, columns 7-8 by two-loop for comparison. 
The size of the required Wilson coefficients c depends on the chosen suppression scale |2]) parametrized by 5, four plausible 
values of which are given for each model: naive or reduced Planck scale, either or not taking the running of Newton's constant 
into account, cf. Table [I] For the last column, see below estimate ( |13[ ). 



These unification results largely agree with [15] . There is, however, some disagreement on a sensible normalization 
among the dimension-5 operators with Higgses in different representations, which directly affects the required sizes 
of the Wilson coefficients c. We have chosen to write all possible dimension-5 operators in a common form with 
uniform normalization conventions (see the Appendix), and believe that this allows for a sensible cross-comparison 
between the Wilson coefficients of different operators. Furthermore, within these conventions, |c| ~ 1 is believed to 
be a natural size for the dimension-5 operators from an effective field theory point of view (see Section [Tl] and the 
Appendix). Since the low-energy gauge couplings a s (mz) had not been measured well at the time and were therefore 
partly treated as free parameters, Table [TTT] cannot be directly compared to the pioneering work that looked 

only at the case of a sole 24 Higgs, but the results tend in the same direction. 



Table III also shows results obtained by using two-loop renormalization group equations [17] (neglecting matter 



fields) instead of (A.30). A full two- loop treatment would further incorporate the effect of thresholds, as partly done 



in |15j , and couplings to the fermion and Higgs fields, increasing arbitrariness and decreasing predictivity |23l | 
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FIG. 1: Numerical gauge coupling unification results for the SU(5) model with a 24 and a 75 Higgs. Left panel: For different 
fixed ratios V24 : «75 of the Higgs vevs, the curves characterize the size of the Wilson coefficients c; necessary to achieve 
unification at Mx- The straight horizontal dashed and solid lines indicate two possible choices for the parameter £ run = 1.6 
and ^'ed = 8.0 in ^ (cf. also Table l|, and the vertical lines the corresponding suppression scales Mpi = 1.2 x 10 19 GeV/£ (it 
may not make sense to consider scales Mx larger than Mpi). Right: The e s in (JsJ) necessary to achieve unification at a given 
Mx (not dependent on the ratio V24 : W75 of Higgs vevs). 



there would also be some renormalization scheme dependence. The one-loop values are not significantly altered, so 
a one-loop analysis can show, with reasonable numerical accuracy, whether the modified condition ([6| can make the 

standard model gauge couplings unify, as a proof of principle. 

Furthermore, endowing the measured low-energy couplings ( |A.32 ) with error bars would put uncertainties on the 
Mx and c necessary to achieve u nificat ion under these initial conditions; for example, in the case of a 75, the error 
bars (less than 4% [Hj) given in \A 32\ result in M x = (8.1 ± 6.9) x 10 15 GeV and c = (-129 ± 93)/£, which would 
not significantly influence whether (natural) unification is considered a reasonable possibility or not. (In a previous 
paper ifv, however, the relative effect from the uncertainty on the necessary Wilson coefficients c was seen to be much 
bigger, as very small Wilson coefficients are required in supersymmetric unification models in the first place.) Besides 
being predictive, our analysis shows the effect due solely to the modified unification condition ^ cleanly separated 
from other effects like, e.g., threshold corrections [231 [2"4"] . 



B. Unification in models with two dimension-5 operators 

In the case of multiple dimension-5 operators ([l} the modification to the gauge kinetic terms Q is a linear combi- 
nation of the effects from single operators, weighted by the Wilson coefficients and Higgs vevs, see ((Hi. The idea is 
that now, as one can vary these contributions with the model parameters in a continuous way, one can find subsets 
of parameter space that result in perfect gauge coupling unification at continuously variable scales Mx ■ 

In the case of two dimension-5 operators (not in the same representation) these subsets are two-dimensional, by 
counting the number of parameters (cj and Vi for each Higgs i) less the two unification constraints (JgJ) : 2-2 — 2 = 2. 
To emphasize the possibility of varying M x continuously, we choose to parametrize these subsets by the unification 
scale Mx itself and by the ratio v\ : V2 of the two Higgs vevs, which, without fine-tuning the Higgs potential, should 
probably be within an order of magnitude of each other. 

Numerical results for the model with a 24 and a 75 Higgs are shown in Fig. [T] For a given choice of maximally 
(naturally) acceptable Wilson coefficients max{|c24|, | C75 1 } , the lowest unification scale M x can be achieved for the 
(not finely-tuned) vev ratio of roughly 1*24 : W75 = 1:3 (solid bold curve in the figure). For example, demanding 
|c24 1 , IC75 j < 1 and choosing the suppression scale Mpi = 1.2 x 10 19 GeV/£ ro d = 2.4 x 10 18 GeV, then any unification 
scale M x > 5 x 10 17 GeV can be achieved, whereas, if one allows | c 2 4 1 , |c 75 | < 5 with M Pl — 1.2 x 10 19 GeV/£™jJ — 



1.5 x 10 18 GeV, then any M x > 3 x 10 16 GeV is possible, getting close to the contraint on proton decay (Section IV I 
The right panel shows the modifications to the gauge kinetic terms Q necessary to achieve unification at Mx in this 
model; in no case are big e s required that would invalidate the analysis (in particular, e s > — 1 always). 

Similar results obtain for the two models with 24, 200 resp. 75, 200 Higgses, see Fig. [2j Within the allowances 
of the previous paragraph, the widest unification scales achievable are Mx > 1.4 X 10 16 GeV, namely for the model 
with 75, 200 Higgses, when accepting max{|c75|, |c2oo|} < 5 and assuming Mpi = 1.2 x 10 19 GeV/^™. Requiring 
maxi \a\ < 1 and assuming the non-reduced Planck scale, one can achieve unification at any Mx > 8 x 10 17 GeV 
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FIG. 2: Similar to the left panel of Fig.fi] but for the models with 24, 200 Higgses (left; £ Iun = 2.42, = 12.1) and 75, 200 
Higgses (right; £ run = 2.68, ££ed ~ 13.6)7 For the 24, 200 model, the ratio of Higgs vevs that achieves any chosen Mx with 
the smallest Wilson coefficients is roughly V24, : «200 = 1 : 3 (solid bold curve), and for the 75, 200 model this ratio depends on 
Mx and lies between 1 : 3 and 3:1, requiring no fine-tuning of the Higgs potential. In both models, \e s \ < 0.5 at all unification 
scales Mx shown here (cf. right panel of Fig. []]. 



(roughly the same for both models). 

In models with two dimension-5 operators ([I]) involving Higgs multiplets in identical irreps (e.g. two 24 Higgses), 
one cannot shift the unification scale continuously since the ratio between the e s cannot be varied continuously. In 
fact, only the unification scales Mx given in Table [ITT] are possible, with Higgs vevs and Wilson coefficients modified 
by O(l) factors. Furthermore, in models with one singlet and one non-singlet Higgs (e.g., one 1 and one 24 Higgs) 
the unification equations do not have physically sensible solutions (formally, the equations yield a s (Mx) ~ at any 
M x ). 



C. Effects in general models 



The behavior from Figs. [T] and [2] is generic for models with at least two dimension-5 operators, as we will see now. 
In particular, there are no points in parameter space which yield unification at some given Mx for significantly smaller 
Wilson coefficients than shown in these figures. This is true in models with at least two dimension-5 operators ([I]) 
and involving Higgses in at least two unequal irreps from the set 24, 75, 200; in all other models, acceptable exact 
unification via mechanism ([6]) is impossible for general Mx, as we have just seen, and can at most happen at the 
discrete scales Mx from Table [Hl| 

Writing c = X)i=i °i an d v = Si=i v i with h the number of Higgs multiplets at the grand unification scale, equations 



^ 1 ® , (|8| and Q together yield for the Wilson coefficients Cj 



12 



\a a (M. 



x) 




(for any s = 1, 2, 3) 



(9) 



(Higgs fields in representations ^ {1,24,75,200} can easily be accommodated by setting Cj 



0). Then 



setting Vi/v w l/h (no large hierarchies between vevs Vi), maxj \ci/c\ > 0(1 — 5)//i > 1/h (the constant of propor- 
tionality depends on the h iera rchy between and the signs of the c,), an (My) > 1/50 (typical for non-supersymmetric 
unification; see also Table "~ 
gives 



IHj) and J2j c 2{rj) = hC 2 > bh (see Table |Vj), and using that ([9j is valid for all s = 1, 2, 3 



max Cj 



> 



> 




(10) 



The last factor depends on the group theory constants 5 ^ associated with the embedding of the standard model into 
SU(5) and numerically equals (10/V168) -1 (see Table |ll[); the other factor largely characterizes how well (or how 
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FIG. 3: For a model with the three (h = 3) Higgs multiplets 24, 75 and 200, a random sample is taken of 30, 000 points {c;, Vi} 
of the parameter space which yield gauge coupling unification. The point plot on the left characterizes the required sizes of 
the Wilson coefficients C; to achieve unification, as in Figs. [I] and [2] unification at smaller scales Mx requires proportionally 
larger Wilson coeffi cien ts. To illustrate the validity of l]13|, the right panel shows, for each sample point, the O(0. 1) constant to 



achieve equality in ( 13 1; numerical gauge coupling unification at any 10 13 GeV < Mx < 10 20 GeV can be achieved by choosing 
0.5 as this constant, and a lower bound is obtained by the choice 0.1. However, this should not be mistaken: unification at 
small M x < 10 16 GeV requires large dimensionless Wilson coefficients Ci (see left panel). 



badly) the actual running couplings a s (/i) of the standard model unify without any modification to the unification 
condition since it can be estimated with ([6]): 



E 



a G (M x ) 



a s (M x) 
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1 M Pl 
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(for 10 15 < Mx/GeV < 10 



19 \ 



(11) 



(12) 



This estimate captures unification according to the modified unification condition (|6]) pretty well for models with 
several Higgs multiplets (as specified above), as can be seen from numerical studies: from the numerical behavior, 
to achieve equality in (12 1, 1/15 should be replaced by some factor of order 0(0.08 — > 4), depending on the Higgs 



content of the model, on the ratio of vevs and Cj's, and on M x . The following lower bound can be strict in some 
models, but has leeway in most situations: 



max Ci 



> 



O(0. 1) Mpi O(0.1) 1.2 x 10 19 GeV 



M 



x 



(for all < Mx/GeV < 10^ u ) 



(13) 



Furthermore, unification at any given M x can be achieved in any such model when allowing coefficients c-i of the size 
(13) with the constant being O(0.5). 

For justification of (13) and its O(0. 1) factor in a model with three Higgs multiplets, see Fig. [3j especially the 



right panel where the OjO.l) factor to achieve equality in (13) has been plotted for 30,000 randomly sampled points 
of the parameter space that yield unification. (One can formally achieve small Wilson coefficients by making the 
number h of Higgs multiplets large, but these models face issues such as Landau poles close to the unification scale 
and small Higgs vevs which lead to scalar masses potentially far below the unification scale.) An estimate like (13) 
allows one to quickly judge whether, for some given Higgs content, unification via (|6| at a desired scale M x is possible 
with natural-size Wilson coefficients Ci. Models containing Higgs multiplets in representations r ^ {24, 75, 200} have 
larger required Cj than models without such multiplets, see ([9]). Estimate (13) is also quite accurate at the unification 
point M x in models with only one dimension-5 operator, see last column of Table |TO} however, in those models, M x 
cannot be continuously shifted to other values as in the models considered in this subsection here. 

The ingredients of the unification mechanism discussed in this paper are well illustrated by ( |13[ ): For models 
with Higgs multiplets in at least two of the irreps 24, 75, 200, unification at continuously variable scale M x may 
be achieved due to the presence of gravitationally induced effective dimension-5 operators in the grand unified 
theory, since two or more such operators allow for a continuous set of solutions to unification condition ([6]). To 
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actually achieve unification at a given Mx, however, these operators need to have the right sizes Cj, see (13) for an 
estimate, and these sizes influence whether one considers unification at Mx possible in a natural way (requiring, very 
roughly, |cj| ~ 0(0.1 — > 10) or so). The sizes of the Cj are further directly affected by the choice of the Planck scale 
Mpi ^ that suppresses the dimension-5 operators 0. Reasonable choices in (13) are £ = 1 (naive Planck scale) 



or, more commonly, £ rcd = \/8tt s» 5 (reduced Planck scale), possibly further enhanced by factors 0(1.5 — » 5) if the 
running of Newton's constant ([3| in such models is taken into account as another, additional gravitational effect (cf. rj 
in Table [I]). These latter considerations about the appropriate choice of the Planck scale merely influence the required 
numerical values of the Wilson coefficients Cj by factors of order 0(1 — > 25), and thus are only secondary to the 
possibility of continuous variability of the unification condition ([6]), which facilitates exact gauge coupling unification 
for some values of the Ci in the first place. 

In this section we have discussed how models can achieve exact gauge coupling unification through unification 
condition ^ and whether this is possible with natural coefficients. The former is a formal numerical requirement 
which grand unified theories ought to obey. But it is a further question whether those models are physically sensible 
and viable, e.g. if they are allowed by experimental constraints. In the following two Sections |IV| and |V} we will look 
at some of the physics of these unification models. 

IV. NON-SUPERSYMMETRIC UNIFICATION NEAR THE PROTON LIFETIME LIMIT 

In this and in the following section we will demonstrate and analyze two possible implications of models which 
employ the unification mechanism just described. 

Here, we display and examine models which allow non-supersymmetric unification of the standard model into SU(5) 
without intermediate symmetry breaking scales, while easily avoiding constraints from proton decay, contrary to a 
widely held belief. This scenario seems attractive as it does not require presently unobserved supersymmetry and 
as non-supersymmetric models offer the possibility of describing physics up to the onset of quantum gravity in a 
perturbative way, see the right column in Table [TJ and the last paragraph of Section [TTj However, the proposed models 
are not completely minimal in their Higgs content as mechanism ([6]) requires at least two multiplets for continuously 
shifting Mx around (the minimal choice being a 24 and a 75 Higgs), whereas one 24 alone would already be sufficient 
to break the grand unified symmetry; nevertheless, as described in Section [TTJ when this Higgs content is assumed, 
the presence of the gravitational operators ([!]) is to be expected and does not further detract from minimality. 

In non-supersymmetric models, proton decay is mediated by gauge d — 6 operators (baryon number violating oper- 
ators after integrating out the superheavy gauge fields) and Higgs d = 6 operators (after integrating out superheavy 
Higgses) . The latter are strongly dependent on the Higgs sector and Yukawa couplings; since they are, apart from the 
potentially dangerous triplet route, generically less important than the gauge d = 6 operators (see [34J for a review) , 
we will concentrate on the gauge contribution for the following estimates (also neglecting potential "textures" in flavor 
space that could partially rotate away the effective gauge d = 6 interactions, yielding a slower decay rate and weaker 
bounds [35]). Under these assumptions, the proton decay rate l/r^^^^o) through the dominant (in our models) 
channel p — > e + ir° is related to the superheavy gauge boson masses Mx ^ by 

M 4 

r (p ^ e+7r o) = 0(1) , (14) 

with the proton mass m p and we take the O(l) constant of proportionality to be 1 for the following. The current 
experimental bound on this decay channel is [20] 13S] 

T( p ^e+K°) > / • 10 33 years with / = 1.6 or / = 5.4 . (15) 

This sets a proton lifetime bound on the unification scale Mx in non-supersymmetric models: 

/ / f \ 1/4 

M x > (40a G ) 1/2 ( 2.4 x 10 15 GeV « (1.9-3.8) x 10 15 GeV , (16) 



where the latter illustrates the range given by (15) and a G = 1/60 — 1/30 (typical for non-supersymmetric models 
using M, cf. also Table |lll|. Within the next ten years, if proton decay remains unobserved, improvements in the 



bound (15) up to / = 100 are expected [37], constraining Mx > 8 x 10 15 GeV. 



One can now see that all of our models that achieve natural gauge coupling unification, i.e. with Wilson coefficients 



of order |ci| max < O(10), satisfy the proton decay constraint (16) easily, and are also fairly safe against the expected 
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ten- year improvements in the bound. For example, of the models with one Higgs multiplet (see Table III), only the 
200 model (and possibly the 75 model, depending on the choice of Mpi) feature natural gauge coupling unification, 
with the unification scale 5 x 10 18 GeV (resp. 8 x 10 15 GeV) above the bound ( 16 ) in either case. This also holds for 



models with two Higgs multiplcts, see Figs. [T] and [2] Under the condition |cj]^ax < 5, the lowest unification scale 
Mx = 1.4 x 10 16 GeV is here achieved in the model with a 75 and a 200 Higgs (Fig. [2J right panel) when assuming 
£ = Ccd = 13-6; the other two models with two multiplets automatically have Mx > 3 x 10 16 GeV if gauge coupling 
unification with Wilson coefficients of this size is required. 

Quite generally, any reasonable non-supersymmetric unification model featuring natural gauge coupling unification 
via ^ satisfies the proton decay constraint (16). This can be seen from the general estimate (13) of the parameters 
necessary for gauge coupling unification: 

M X > °M 1 - 2X1 ° 1 , 9GCV > -^ L2X1 ° 19GCV = 4xlO-GeV, (17) 
£y/h maxilcil ~ 15^4 10 



so that ( 16 1 still holds despite all the very conservative parameter choices in ( |17| which are by no means necessary 
or particularly desirable. Models with more natural parameters (i.e., with smaller |cj| max , reasonable h, and possibly 
choosing £ re( j = v°7r) satisfy the proton lifetime constraint very easily. 

The intuitive reason why our models naturally obey the proton decay constraint so readily is clear: In non- 
supersymmetric models, the couplings a s (/j,) miss each other by quite a bit, so the modifications to the unification 
condition ^ need to be relatively sizeable e s ~ O(0. 1) in order to achieve exact gauge coupling unification; then, 
since we require naturalness |ci| max < 0(10), these sizes e s ~ CiMx/Mpi of the effective gravitational corrections Q 
must be due mainly to the proximity of the Planck scale Mpi = 1.2 x 10 19 GeV/£ to the unification scale Mx and, 
consequently, to the gauge boson masses ((8|). Typically, such heavy gauge bosons ensure the proton decay constraint 



{16). 

On the other hand, whereas the natural models easily satisfy current proton decay limits, there are situations 
conceivable in which the constraints in the not too distant future will come close to testing some of the models, either 
excluding them or, if the proton is actually seen to decay, strongly restricting their parameter space. For example, a 
model with Higgs content 24, 75 and 200 can achieve unification at Mx = 8 x 10 15 GeV, the projected constraint 
in less than ten years from now [37], if merely the reduced Planck scale (£ = £ re d = v87r ps 5) and somewhat large 
coefficients |ci| max ~ 20 are accepted (see Fig. [3]). Alternatively, the smallest model with two Higgs multiplets 24 and 
75 can reach this limit when taking £ = = 8 and allowing |ci| max ~ 15, and similarly can the other models with 
two multiplets (in each example, only the product £|ci| max is fixed whereas the sizes of both factors can be traded 
back and forth). Regarding the exact numerical values in this discussion, one has to keep in mind that there is some 



uncertainty in the exact numerical O(l) factor in (14), although it enters the final bound (16) only to the power of 
1/4, and in higher-order and threshold corrections to the gauge coupling running. 

We have seen that the natural models with gauge coupling unification through ^ are physically viable as they 
automatically evade the proton lifetime constraint that naively excludes non-supersymmetric unification. However, 
some of the less natural unification models will come close to the proton limit in the near future, which will exclude 
some of them, or constrain their parameter space drastically if proton decay is observed. 



SIMULTANEOUS UNIFICATION OF GAUGE AND GRAVITATIONAL INTERACTIONS 



As can be seen from Figs. [IH3]or from estimate (13), gauge coupling unification at large Mx ~ Mpi is possible 
through mechanism ^ naturally, i.e. with small Wilson coefficients Cj. Unification of the three standard model gauge 
interactions with each other on the one hand and with the gravitational interaction on the other hand at the same 
scale ("gauge-gravity unification"; cf. [5TJ [55]) therefore constitutes another scenario naturally achievable through the 
modified unification condition ^ . However, one should keep in mind that our analysis is based on an effective theory 
approach and we should therefore consider carefully the expansion we are using. If the operators of dimension 5 and 
higher that we are considering are of non-perturbative nature, then the expansion is in powers of l/Mpi and one may 
worry that it could break down once the Higgs vevs get close to the Planck mass Mp[. In any case, it is an interesting 
numerical coincidence that the corrections to the unification condition allow to shift the unification scale close to the 
Planck scale for a natural set of parameters. 



Approximate parameter values Cj of such gauge-gravity unification models can be read off from estimate ( 13 ) by 
setting M x = M PI /0{1): 

max| Ci | « 0(1) Q( °- 2 1] , (18) 
* vh 
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where the 0(0.2 — > 1) estimate stems from the Mx ~ Mpi range of numerical studies like Fig. [3] (right panel) and h 
denotes the number of Higgs multiplets at the grand unification scale. This suggests the possibility of very natural 
gauge coupling unification at or near Mpi and is already apparent in the only model with one Higgs multiplet that 
achieves unification close to the Planck scale: The 200 model has Mx = 5.2 x 10 18 GeV, which is related to the 



Planck scale by Mx — 0.4£Mpj, and the required Wilson coefficient is natural c = 0.53/£ (see Table III I. As another 
example .gauge coupling unification can be achieved for appropriate parameter choices in the model with a 24 and a 
75 (Fig.llj at the scale Mx = Mpi = 1.2 x 10 19 GeV/£ for any of the reasonable exemplary choices £ = 1, £ ro d = V^tt, 
^run _ ^ q Qr £run _ g q. correS p nding model parameters |ci| max = 0.23, 0.21, 0.23 and 0.20 are all natural and 



of almost equal sizes independent of £, see (18). Similar numerical estimates hold for the other two- and three-Higgs 
models. In these models, any O(l) factor in the relation Mx = MpijOiV) can be easily accommodated as well by 
corresponding O(l) changes to the Wilson coefficients, whereas this is not possible for the sole-200 model as it only 



permits one discrete unification scale (Table III I . 

In the effective field theory spirit of SectionjllJ operators of dimension higher than 5 are also present, e.g. higher- 
dimensional generalizations of |l]): 

C = -^tHiHzG^ + ^ Hl H 2 H 3 G^G^ + ... . (19) 

After the Higgs multiplets acquire vevs at the scale Mx = Mpi/0(l), they contribute to the gauge kinetic terms Q 
as well: 

1 - j (i + u + ^ + w + ■ ■ ■) n^r ■ » 

with the corrections e s ~ c^SgMx/goMpi from the dimension-5 operators ([5]). Depending on the group theory factors 
df'' 7 ^ (analogous to the S a in Table [ll| and on the constant in the relation Mx — MpifO(l), this expansion might 
or might not be controlled perturbatively. If it is not, one cannot claim perturbative gauge-gravity unification at the 
Planck scale. Nevertheless, the fact that mechanism Q allows in principle to naturally adjust the unification scale 
to a high scale ~ Mpi might at least be taken as a hint that gauge-gravity unification is a possible scenario, even if 
the necessary parameter values or the last piece of the evolution cannot be computed perturbatively. 



VI. THE NON-SUPERSYMMETRIC SO(10) CASE 

In this section we give an account of the effects from the dimension-5 operators in models with grand unified 
group G — 5O(10). We will describe the setup for 5O(10) by emphasizing the differences to the SU(5) case described 
in Section [TT] and give a few numerical unification results similar to Section |III| We find that the effects from Sections 
IV (on proton decay) and [V] (on gauge-gravity unification) can occur with similar sizes for 5O(10) as well, although 



the 5O(10) formalism is more general than the 5C/(5) one; in particular, a continuously variable unification scale 
Mx can now be achieved with only one single dimension-5 operator. For this reason, and for the beauty of 5O(10) 
unification, we find this treatment worthwhile. 

For G — 5O(10), there are two inequivalent ways to embed the standard model group G321 = SU{2i)c x SU(2)l x 
U(1)y into G consistent with the charge assignments of the standard model fermions: the "normal embedding" 
G321 C SU(5) C 5O(10) [35], and the "flipped embedding" G 321 C 517(5) x U(l) x C 5O(10) with G 32 i £ SU(5) 
[35"] (see the Appendix for a careful treatment of the 5O(10) group theory). Statements without qualifier in the 
following apply to either embedding. Note, nowhere are we implying that either SU(5) or SU(5) x U(l)x (or any 
other subgroup) be intermediate unbroken symmetries at any scale. 

In the 5O(10) grand unified theory, the gravitationally induced effective dimension-5 operators ([I]) can be formed 
with Higgs multiplets Hi in any of the four irreps ri = 1, 54, 210, 770. However, contrary to the SU(5) case, the 
vev (Hi) is not necessarily uniquely specified (up to normalization) by requiring it to be invariant under the standard 
model group G321; rather, this G32i-invariance merely restricts to a 3- resp. 4-dimensional subspace inside the 210 
resp. 770 irreps, so that the length as well as the direction of the vev has to be specified for Higgses in these irreps. 
This will have the important consequence that with merely one dimension-5 operator ([T]), built with either a 210 or 
770 Higgs, one can continuously vary the ratio among the e s in the unification condition Q by continuous variation 
of the vev direction, ultimately leading to a continuously variable unification scale Mx', the SU(5) case could achieve 
this continuity only with at least two Higgs multiplets. 

To parametrize these vev directions for calculations, one has to specify a basis in these 3- or 4-dimensional subspaces. 
Two possible, distinct choices of this basis are given in the Appendix: the vectors in the first resp. second basis have 
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definite transformation properties under the SU(b) x U(l)x resp. SU(A)c X SU(2)l x SU(2)r maximal subgroups 
of 50(10) (see the first two columns of Tables [VlT resp. |VIII| for these transformation properties). We choose to 
classify according to these two subgroups since they are the only maximal subgroups that can occur as intermediate 
symmetries in 50(10) — > G321 breaking consistent with the standard model charge assignments. These two choices 
therefore readily facilitate the analysis of 50(10) breaking with intermediate gauge symmetries, see below (the absence 
of intermediate symmetries is assumed for now). 

With such a chosen basis, the modifications e s to the three standard model gauge kinetic terms Q are now, instead 
of j5j, double sums ( |AT28| > 



e * = E^£w^ (for , = 3,2,1), 



(21) 



where the second sum runs over the directions j = 1, . . . , 3 (4) if Hi is in the 210 (770) representation (for the 1 and 

54 irreps, j = 1 only). The Clebsch-Gordan coefficients 5* depend on the embedding G321 C 50(10) (normal or 
flipped) , on the representation of the Higgs Hi and on the choice of basis vect ors j in the standard model singlet 
subspace of the respective representation. These 5« can all be taken from Table VII (for basis vectors with definite 
5 £7(5) x U(l)x transformation properties) and Table VIII (definite SU(4)c x SU(2)lX SU(2)r transformations) via 



One other difference between the 50(10) and the SU(5) cases is the fact that none of the 54, 210 or 770 multiplets 
can give mass to all of the non-standard model ("superheavy") gauge bosons and that, furthermore, the masses of 
the gauge bosons which actually do acquire mass from these multiplets are generally unequal. Therefore, both the 
statement from the SU(5) case that the mass M g b of the superheavy gauge bosons may originate solely from Higgses 
involved in dimension-5 operators ([I]) and the requirement ^ that this mass be close (or equal) to the unification 
scale Mx have to be amended. First, additional Higgs multiplets in other irreps (typically 45, 16, or 126; see also 
Table [l| are necessary to break the 50(10) gauge symmetry down to the standard model G321 and they contribute 
to the masses of the superheavy gauge bosons as well. And for definitcness, we now assume that the Higgs multiplets 
involved in dimension-5 operators ([!]) account for half of the average squared superheavy gauge boson masses, and 
we require that this averaged superheavy mass now be equal to the unification scale, see (A. 25 1 in the Appendix. 



The necessity for some such choice causes uncertainty in actual numerical unification calculations greater than in the 
SU(5) case. 

As mentioned above, continuously variable unification scales Mx can be achieved in 50(10) models with a single 



Higgs in cither the 210 or the 770 representation, since the direction {vu = uj}j— 



continuously in (21) 
attractive 



(*=l)j/i=l,...,3(4) 



of their vcv can be varied 



Minimizing the Higgs content and the number of dimension-5 operators in this way seems 



Further note, that a 770 causes the strong coupling regime of the grand unified theory to be quite close 
to the unification scale (Table [T]), although, depending on the gap between the unification and the Planck scale, the 
theory may still be perturbative up to the onset of quantum gravity. 

As a numerical example, Fig. [4] shows gauge coupling unification results in the 50(10) model with a single 210 
Higgs (identical results obtain for both the normal and the flipped embedding due to the orthogonal relation ( |A.18 )), 
similar to Fig. [3] for SU(5) with three Higgses. The lower bound (131 holds here as well, even with the same factor 
O(O.l), and unification can be achieved for any Mx in the displayed range if this factor is allowed to be 0.5 (see right 
panel). Therefore, unification is naturally safe from the proton lifetime limit in 50(10) as well, but can also come 
close to it for some points in the parameter space, and unification near or at the Planck scale is achievable in a very 
natural way, both similar to the corresponding 5t/(5) scenarios described in detail in Sections IV and|v} 



In a scenario with several breaking scales, if the intermediate gauge symmetry is 5t/(5) x U(l)x (or merely SU(S)) 
or the Pati-Salam group SU{A)c x SU(2)l x SU(2)r, respectively, the directions of both the 210 and 770 vevs at 
the unification scale Mx are fixed, as th ese vevs have to be singlets under the intermediate group, which allows only 

Therefore, similar to the case of an SU(5) 



the 5s from rows 1 



3, 6 in Table 



VII 



VIII 



resp. rows 1, 2, 3, 6 in Table 

unified group, only linear combinations of the effects from at least two operators ([!]) can achieve continuously varying 
Mx ■ The effect of a dimension-5 operator with a Higgs that acquires its vev at an intermediate scale Mj < Mx is 
suppressed relatively by ~ Mj/Mx- Note, since there are no singlet vevs for a 54 Higgs under the SU(5) intermediate 
group, it can acquire nonzero vev only at the final breaking scale to the standard model group G321 and is therefore 
almost irrelevant for our effect in the case of an SU (5) intermediate gauge symmetry. 

In several previous analyses of the 50(10) case (e.g., [H]), only Higgs vevs that are invariant under intermediate 
gauge groups larger than G321 have been considered (namely under the SU(5) Georgi-Glashow or the SU(4)c x 
SU(2)l x SU(2)r Pati-Salam subgroups), thereby omitting most of the possibilities discussed above for 50(10) 
breaking without intermediate symmetries, i.e. the continuous variability of Mx with merely one dimension-5 operator. 
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FIG. 4: For the 5*0(10) model with a 210 Higgs, a random sample is taken of 30,000 points {C210, f(2io)j} °f the parameter 
space which yield gauge coupling unification, cf. Fig. [3] Estimate ( 13 1 is valid here as well (see right pan el) , and numerical 
gauge coupling unification at any 10 13 GeV < Mx < 10 GeV can be achieved with the O(0. 1) constant in ( 13 1 being at most 
0.5 (the caveat from Fig. [3] applies here as well). 



Exactly analogous effects obtain for other grand unified gauge groups, like -Eg, and they are even of similar size as 
in SU(5) and 5O(10) [HI [25]. Furthermore, the formalism necessary for the 50(10) treatment is already the most 
general for any unification group, so we leave the analysis of gauge coupling unification through the mechanism of 
dimension-5 operators (JT1) at the S77(5) and 5O(10) cases discussed so far. 



VII. COMPARISON TO THE SUPERSYMMETRIC CASE 



The influence of the dimension-5 operators has been studied 2 in models of supersymmetric grand unification 
[T6T - H8] as well, although most studies have focused on the effect of the supersymmetrized version of in the creation 
of non-universal gaugino masses in a scenario where the UV completion is a M = 1 supcrgravity theory, see [25] and 
references therein. The group theory needed to calculate the resulting gaugino mass ratios is very similar to the 
formalism for obtaining the modified unification condition ^ [5SJ[27]; see the Appendix for the detailed SU(5) and 
5O(10) group theory. 

The effect of these dimension-5 operators on the unification of gauge couplings in supersymmetric theories and on 
the unification scale has been noted before [2] [5] [TH] HH [37]. In some of these works, however, the low-energy 
inputs a s (mz) or the supersymmetry breaking scale msuSY (above which the /3-function coefficients (A. 31 1 are taken 
to be (pi, hi, 63) = (33/5, 1, —3)) have been treated as uncertain parameters, and mostly only the effect from a single 
dimension-5 operator had been taken into accou nt (exc ept in [T3"]). 

Fixing the low-energy gauge coupling values ( A.32[ ) and msuSY = 1 TeV [T5], as roughly required to solve the 



hierarchy problem of the standard model, one obtains for supersymmetric £77(5) models with one dimension-5 operator 
the unification possibilities in Table IV (at one loop), in analogy with Table [TTT] for non-supersymmetric SU(5) (note, 
supersymmetric SU(5) with a 200 Higgs enters strong coupling shortly above the unification scale, cf. Table [I]). 
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1/qg 
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max s | e s \ 
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24 


1.1 x 10 16 


25.9 


31.1/f 


2.4 x 10 16 


0.024 


1 


5 


2.3 


11.6 


75 


3.2 x 10 16 


25.8 


-12.0/C 


5.6 x 10 16 


0.033 


1 


5 


3.1 


15.4 


200 


8.5 x 10 16 


26.8 


13.1/* 


1.2 x 10 17 


0.105 


1 


5 


4.4 


22.1 



TABLE IV: This table shows the parameters that a SUSY-SU (5) model must have if, at one loop, unification happens by 
means of mechanism Q with only one dimension-5 operator |l]). The sizes of the required Wilson coefficients c depend on the 
chosen suppression scale parametrized by £. 

For a general number of Higgs fields, the situation is shown in Fig. [5] in analogy to Fig. [3] As can be seen, 
exact supersymmetric gauge coupling unification around 2 x 10 16 GeV can be achieved with natural-sized dimension-5 
operators ([T]), which is expected as this is where the standard model gauge couplings come closest to each other in 
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2 06 




5: Similar to Fig. [3] but for the supersymmetric SU(5) model with the three (ft, = 3) Higgs multiplets 24, 75 and 
Again, numerical gauge coupling unification at any 10 13 GeV < Mx < W 20 GeV can be achieved by choosing 0.5 as 



FIG. 
200. 

the constant in (131 (but cf. caveat from Fig. [3| . Most of the randomly chosen unification points lead to unification scales 
Mx ~2x 10 16 GeV, where, in supersymmetric models, the standard model gauge couplings almost meet. 



supersymmetric models; but higher unification scales are also possible in a natural way with certain values of the 
coefficients Cj. Note that the supersymmetric unification model from Fig. [5] has a Landau pole roughly half an order 
of magnitude above the unification scale Mx , cf. Table [TJ 



VIII. SUMMARY AND CONCLUSION 



We have studied the effects from multiple gravitationally induced dimension-5 interactions c i H i G lljV G llv /AM P i (HI), 
naturally present in any grand unified theory, on the unification of gauge couplings. These operators can mod- 
ify the canonical gauge kinetic terms Q and the coupling unification condition (JsJ by up to e s ~ Ci{Hi)/Mpi ~ 
CiMx/gaMpi ~ 0.1% - 30% (5) after the Higgs multiplets acquire vevs in grand unified symmetry breaking. The 
size of the effect depends on the unification scale Mx, the Wilson coefficients 0.1 < |e,j| < 10, and the suppression 
scale 5 x 10 17 GeV < M Pl < 1.2 x 10 19 GeV @. 

Modifications e s of this size can effect perfect gauge coupling unification, at scales Mx significantly different than 
naively expected, as the standard model gauge couplings a s (/x) (s — 1, 2, 3) differ from each other by at most < 50% 
in the wide range 10 13 GeV < fj, < 10 19 GeV, both in the case with and without supersymmetry. When two or more 
dimension-5 operators are present in a theory (or a single 210 or 770 in 50(10), see Section VI), then the scale Mx 
of gauge coupling unification can be varied in a continuous and controlled way as a function of the Wilson coefficients 
d and Higgs vevs (Hi), which are practically model parameters. We have focused on non-supersymmetric 5t/(5) and 
5O(10) theories, for which grand unification seemed, previously, to be difficult [12], and we find that gauge coupling 
unification is possible, in a natural way, at any Mx ~ 10 17 GeV — 1.2 x 10 19 GeV. For illustration we now give a few 
numerical unification results (all non-supersymmetric models): 

(a) SU(b) with 24 and 75 Higgses, see Fig. [lj Perfect gauge coupling unification at Mx — 10 17 GeV happens 
for some choice of parameter values C24, C75 with \a\ < 3, if Mpi = 1.2 x 10 19 GeV/y/8ir = 2.4 x 10 18 GeV is 
assumed; the modifications e s to the gauge kinetic terms are all less than 15%. 



(b) 5t/(5) with 24, 75 and 200, Fig.|3]and estimate (T3}: Unification at any M x > 10 18 GeV is possible for several 
choices of C24, C75, C200 with |c;| < 2, when Mpi = 1.2 x 10 19 GeV is assumed. In fact, the bigger a unification 
scale one wants to achieve, the smaller Wilson coefficients are required: maxj |cf| rs 2- (10 18 GeV /Mx); analogous 
statements hold (at fixed Mpi) in all examples (a) - (c) with continuous variability of Mx, cf. (13). 



(c) 5O(10) with a single 210, Fig. [4] Exact unification at Mx = 3 x 10 17 GeV is possible with c 2 io ~ 2.5 for a 
certain (continuous) set of directions of the Higgs vev (H ab ), if M Pl = 1.2 x 10 19 GeV/V&K = 2.4 x 10 18 GeV is 
assumed. 



(d) SU(5) with 24 and 75, Fig. [T and equation (18 1: Numerical gauge coupling unification at the very high scale 
Mx = Mpi, i.e., at the Planck scale, is possible for some choice of Wilson coefficients with | C24 1 , | C75 1 < 0.25 
(irrespective of the specific choice of M Pt ) . 
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(e) SU(5) with a single 200, Table III Exact gauge-coupling unification can happen, but only for a discrete choice 
C200 = 0.5 [15] (here, Mpi — 1.2 x 10 19 GeV is assumed); in this case, the unification scale Mx = 5.2 x 10 18 GeV 
cannot be shifted continuously. 



We have shown (Section IV I that non-supersymmetric unification, achieved via these dimension-5 operators, is not 



in conflict with the current bound on the proton lifetime ( 16 ), since, with natural-sized Wilson coefficients Cj, only 



unification scales Mx > 10 16 GeV are reasonably possible, cf. (17 1. Consequently, experimental improvements of the 



proton decay bound will constrain our models only weakly within the foreseeable future. Our main conclusion is that 
fairly minimal models of non-supersymmetric SU{5) and 5*0(10) unification are easily viable through the mechanism 
described. In particular, supersymmetry does not have to be invoked to save the idea of grand unification. 

We also note that, with two or more dimension-5 operators in the theory, gauge coupling unification at or near the 
Planck scale Mx ~ Mpi can happen for very natural parameter choices Cj ~ O(0.1 — > 1), cf. (18) in Section [vj Of 



course, our approach cannot be fully justified very near the Planck scale, due to incalculable corrections from quantum 
gravity. But the fact that one can, by continuous variation of the Wilson coefficients c,; in a natural domain, push the 
unification scale Mx close towards the Planck scale, may be suggestive of a scenario of simultaneous gauge-gravity 
unification. 

As we have seen, when unification happens according to the mechanism presented in this paper, then unification 
scales Mx much larger than usually assumed are favored, in the sense that probably only those are achievable in a 
natural way. Thus, the mechanism here is distinct from some other known unification scenarios [231 124] . The difference 
lies in the fact that the importance of the dimension-5 operators ([!]) increases for unification scales Mx closer to the 
Planck scale (implying larger Higgs vevs (H)), which causes a bigger modification to the unification condition ^ 
and makes this larger Mx consistent. But, for example, two-loop and threshold effects are not significantly enhanced 
at larger unification scales, since they originate from within the grand unified theory, rather than from intrinsically 
shorter distance effects like strong quantum gravity. 
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Appendix: Normalization conventions and group theory 



In this appendix, we give a self-contained account of the group theoretical aspects and normalization details that 
were omitted in the main text. 

Our first goal is to establish conventions for the quantities appearing in the dimension-5 operators ([I]) of interest, 
aiming towards criteria for assessing the individual size and importance of each operator i based on its Wilson 
coefficient a, and also the relative sizes of different operators ([lj appearing in the same Lagrangian. We intend to 
establish normalization conventions in such a way that operators with Wilson coefficients \ci\ ~ 1 are expected in the 
effective Lagrangian of a grand unified theory after integrating out gravitational interactions, which then helps judge 
the naturalness of a particular model, i.e. its resonableness: If |cj| 3> 1 (for any i), the effective Lagrangian might not 
be a good low-energy description and one might lose perturbative control, whereas |cj| -C 1 constitutes fine-tuning, 
as setting c = does not enhance the symmetry [3] . Note, the calculations and numerical gauge coupling unification 
results in this paper are valid, at face value, regardless of such criteria. 

The effective operators ~ HG^G^" we are considering are formed from two gauge field strengths G MI ,, transforming 
in the adjoint representation G of the grand unified gauge group G, and one Higgs multiplet H , transforming in an 
irreducible representation (irrep). The contraction HG ilv G^ v can yield a gauge singlet only if H transforms in an 
irrep contained in the (conjugate of the) symmetric product (G ® G) s of two adjoint representations of the gauge 
group. For SU(5) and 5*0(10) unified gauge groups G, these possible irreps are the direct summands in the following 
decompositions (for this and other group theory facts, see [2"9"j): 

G = SU{5), G = 24 : (24 fg) 24) s = 1 © 24 © 75 © 200 ; (A.l) 
G = 50(10), G = 45 : (45 <g) 45) s = 1 © 54 © 210 © 770 . 



(In this paper, the 210 of 50(10) will always denote the 210-dimensional irrep in the same congruency class as the 54 
and 770, namely in the congruency class 0, to which also the adjoint G = 45 belongs.) The irreps on the right-hand 
sides of (A.l I are all real representations, as the adjoint G is real itself. 
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Since they are contained in (G <Ei G) s , each multiplet H in an irrep ( A.l ) can be written in component form as H ab 
with symmetric indices a, b of the respective adjoint representation G. Under an infinitesimal gauge transformation 
(1 + ia c T c ) these components then transform according to 



H ab jjab + ia c I yzjaa $bb + $ aa (t G ) bb ) H a b = H ab + ia C [t C G , H]*" , (A.2) 

with the representation matrices (t G ) ab = —if abc of the adjoint G, where f abc are completely antisymmetric structure 
constants of the gauge group. We normalize the real fields H ab and the gauge fields G^„ = d^A^ — d^A® +gaf abc A b l Al,, 
along with the gauge coupling go and the structure constants f abc , in such a way that their kinetic terms have the 
(standard) form 

C = l(d^H ab -ig G A^[t G ,H] ab )(d^H ab -ig G A^[t d G ,H] ab ) - \g%G^ (A.3) 

with the index of the adjoint representation C(G)<5 ab = to{t G t G ) = f acd f bcd normalized to 

517(5) : C(G) = C(24) = 5 , (A.4) 
50(10) : C(G) = C(45) = 8 . (A.5) 



The assignment (A.5 1 follows from (A.4) if one demands the structure constants of 50(10), restricted to its "normal" 



SU(5) subgroup which we will be interested in later, to coincide with the SU(5) structure constants normalized 



according to (A.4), which itself follows from the standard assignment (7(5) — 1/2 for the fundamental representation 
5 of SU(S). We use this standard choice (7(N) — 1/2 for the fundamental representation N of any unitary group 
SU(N). 

In this notation, the unique gauge (and Lorentz) singlet contraction is 

HGG ~ H ab G a ^G b " v . (A.6) 



Writing the dimension-5 operators of interest in this common form for all of the admissible Higgs irreps (A.l), together 
with the normalization ( |A.3 ) of the kinetic terms of the gauge and Higgs fields, establishes a convention (up to the 



sign) for the - a priori arbitrary - overall factor of the singlet contraction of one Higgs multiplet with two gauge fields. 
It moreover asserts, as a convention, that the two operators HiGG and HjGG, even with Higgses Hi, Hj in different 
irreps, are of equal size and strength. The overall sign of the gauge singlet contraction is arbitrary, as it contains 
the Higgs field H to one power and as there is no invariant way to fix the sign of H through the (quadratic) kinetic 
term ( |A.3[ ) . Consequently, this leaves the overall sign of the Wilson coefficients c, in open, but allows naturalness 
assessments based on the absolute values |cj|. 
In we wrote the dimension-5 operators as 



C = y -^Hf -G a uv G b ^ (A.7) 
^ Mp, 1 4 ^ v ' 



(the Hi do not all have to be in distinct representations). The justification for the factor 1/4 is the same as in the 
gauge boson kinetic term (A.3), as the squared time-derivatives of the spatial gauge fields appear as 2(doAf) 2 in the 
contraction -G°„G a » u , 



rereas squared time-derivatives of boson fields canonically appear with factors of 1/2 as 



kinetic terms in the Lagrangian C. To avoid the same overcounting, we explicitly put 1/4 into the operator (A.7) 
as well, noting that this might be a conservative choice (i.e., require a larger value of |cj| in order to achieve the 
same effect). Mpi, which we parametrize as Mpi = 1.2 x 10 19 GeV/£ (cf. Q), should be the scale of the physics 



generating the operators (A.7), and we imagine this to be gravitational interactions (see the beginning of Section 
|ll| ). This completes the setup of normalization conventions that satisfy the goal formulated at the beginning of this 
appendix. 



After a Higgs field Hi has acquired a (nonzero) vacuum expectation value (vev) (Hi), the effective interaction (A.7) 
contributes a term 

/ Tjab\ 

C = A 1 ' G a uu G b ^ (A.8) 
AMpi ^ K ' 



to the Lagrangian, similar to the gauge boson kinetic term in (A.3). Our next goal is to look quantitatively at this 
contribution and its immediate effects. 
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The contributions (A.8) scale roughly like (H°- b ) / Mpi ~ M/Mpi if Hi assumes its vev at energy scale M, and so are 
completely negligible for M anywhere near or below the electroweak scale. On the other hand, above the electroweak 
scale the standard model subgroup G 32 i = SU(3)c x SU(2) L x U(l)y of the grand unified group G is unbroken, so 



the vevs of all scalar fields in the theory, in particular all vevs appearing in (A.8), hav e to be invariant under G321 
(this includes zero vev). This requirement on (77? } constrains the contributions (A.8 1 to the gauge kinetic terms. 
Next we will find, for the cases of 5/7(5) and 50(10) grand unified gauge groups G and for each of the allowed irreps 
(A.l), all possible (Hf) that are invariant under the standard model subgroup G321 C G. 



The method employed is similar to |25| (section I), but we do not only want to find relative contributions, but 
rather also absolute values within our normalization conventions above, and we furthermore want to establish a 
relation between the Higgs vevs and the masses of the superheavy gauge bosons (see later). Let (f> = (4> ab ) be a v ector , 



transforming like (H ab ) as the symmetric product of two adjoint representations G of the gauge group G, see (A. 2 1 
Choose 1 an explicit basis {i^la of generators of the adjoint representation G such that they satisfy the normalization 



(A. 4 ), (A. 5 ) and that the generators of the standard model subgroup G123 correspond to a = 1, . . . , 8 (for the 5/7(3)c 
factor; called "set 7" of the generators), a — 9, 10, 11 (5/7(2)^; set II) and a = 12 (U(l)y; set III or set III', see 
later); this allows for easy examination of the transformation properties of <f> under the standard model subgroup 
G321. Also, the remaining generators a — 13, . . . , d(G) may be identified according to their transformation properties 
under subgroups of the full gauge group G, and we need such a naming scheme later to address the components of 
(f> ab of different superheavy gauge bosons. For our purposes we choose one of the following classifications of generators 
t a (applicable to any representation): 

(a) For 5/7(5): The adjoint G = 24 of SU(5), under which the generators transform, branches under the standard 
model subgroup G321 into: 

24^ (8,1, 0)ffi(l, 3,0)0(1, 1, 0)© [ (3,2,- 1)© h.c. J . (A.9) 



/ 11 in 



IV 



Bcsides the irreps I, II and III, there is only one real representation of G321 contained in the 24, and those 12 
generators make up set IV of generators. 

(b) For 5*0(10): Two different classifications are possible since G321 can be embedded in two different ways (normal 
[55] or "flipped" embedding) into 50(10) with the correct charge assignments for each standard model family 
of fermions from a 16 of 5*0(10). For either embedding, SU(5) in the following denotes the 5/7(5) C 5O(10) 
subgroup which contains the SU(3)c and 5/7(2)/, factors of G321, and U(l)x the Abelian factor such that 
5/7(5) x U{l)x C 5O(10) is a maximal subgroup. 

(i) normal (Georgi-Glashow-like) embedding G321 C 5/7(5) C 5O(10) [35]: 

45 5 ^l ) 24 © 1 © (lOffih.c.) 

°n (8,1,0)0(1, 3,0)©(1, 1,0)© f(3,2,-^)©h. 



1 11 in 




IV 



(3, 2, -) © h.c. J © ( (3, 1, -~) © h.c. ) © ((1, 1, 1) © h.c.) . (A.10) 



6' J V 3 



VIII 



VI VII 



The number of generators in sets IV, V, VI, VII, VIII is 12, 1, 12, 6, 2. Sets I, II, III, IV correspond 
to I, II, III, IV in ( A.9[ ) (this is the Georgi-Glashow embedding). Note, the generators 77/ (generating 



U(l)y C G321) and V (generating U(l)x) both have the same standard model quantum numbers (both 
are standard model singlets), whereas all other sets of generators have distinct ones. 



One way to choose such an explicit basis of adjoint generators is to choose an explicit basis tp for the fundamental representation F of 
G, normalized to tr(ipfp) = C(F)S ab , such that the already obey the desired classification, which may be more easily accomplished 
than the same classification for the t^,. Here, F = 5 or 10 for G = SU(5) or SO(10), respectively, and C(5gu^) = 1/2 and 
C(10 SO (10)) = 1- Then ( t h) bC - -if abc with f abc = -itr ([^,4] *p) /C(F) have the desired properties. 
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(ii) nipped embedding G321 C 5C/(5) x U(l) x C 50(10) with G321 £ 5t/(5) (more precisely, the U(1) Y factor 
of G321 is not contained in SU(5)) |39) : 



45 



S(7(5)xC/(l)x 



24(0) © 1(0) 



ioH/-)eh.c. 



(8, 1,0) ©(1,3, 0)9(1,1,0)® (3,2 



11 



in 



(3,2,--)©h.c. 




VI 



VII 



((l,l,-l)©h.c.) . (A.ll) 



VIII 



Again, only the generators III' and V have the same transformation properties under the standard model. 
The s ets I, 77, IV, VI, VII and VIII of generators here can be identified with the respective sets in 
(A. 10), whereas the two generators t 11 ( generating U(l)y ^ SU(5)) and t v here are linear combinations 
of the generators t 111 and t v from (A.10), 



fin' 



fin 



= V 



with V — 



.lglllglll' 2n/6 gVgUl' 
5 



2^ glll'gV' 



hs v s v ' 



(A.12) 



where s III ,s v ,s 111 ,s v = ±1 are signs that depend on the choice of signs of the generators t 111 , t , 
t 111 and t , which are not fixed by our normalization conventions above. V is a real orthogonal matrix, 



= I 2 , and the gauge bosons F^ 1 , F^ u corresponding to the generators t 



VV A 

gauge bosons F** 1 , F% v via (F IH ',F V ') = (F IH ,F V )V T . 



111 , t v arc related to the 



In this explicit basis, the requirement that cf> be invariant under G321 then, according to (A.2|, translates into 



(t 



GOG 



\(ab)(a'b') ^a'b' 



for all c = 1, 



12 



(A.13) 



Solving this system of (linear) equations (along with enforcing symmetry <\> ah — <f> ba ) yields 4 (resp. 9) linearly 
independent solutions (t = 1, . . . , 4 resp. t = 1, . . . , 9) in the case of G = SU(5) (resp. G = 5O(10)), their 
linear combinations exhausting all standard model singlet vevs contained in (G <g> G) s . In order to find all possible 
standard model singlet vevs for each of the Higgs irreps ( |A.l ) , one has to form linear combinations of the solutions 
4>t that transform in those irreps. To this end, one evaluates the quadratic Casimir operator on each of the linearly 
independent solutions <pt, expressing the result again as their linear combination, 



d(G) 



Co 



[*G.[*G,&]] 



, C ut <t>u 



(A.14) 



which is possible since the restriction of the action of the Casimir operator to an irrep is proportional to the identity 
on that irrep; Table [V] gives the respective constants of proportionality ("quadratic Casimir invariants") for irreps 

with linearly 



of several groups of interest here and later. The matrix (c„t) is diagonalizable, J2t CutV t 
independent eigenvectors v {w \ w = 1, . . . , 4 (9) for G = SU{5) (G = 5O(10)). The 4 (resp. 9) vectors $„, = £ 
then form a basis of the standard model singlet subspace as well, and each <fr w transforms in a representation r of G 
whose quadratic Casimir invariant G 2 (r) is equal to the eigenvalue ■ Since, in the 5C/(5) as well as the 5O(10) 



t "t 



case, the different irreps in ( A.l I have distinct quadratic Casimir invariants, and since each irrep only occurs once in 
the direct sum decomposition (A.l I of (G ® G) s , it is ensured that each <& w actually is contained in an irrep r, and 



r can be read off from the corresponding eigenvalue via Table |_ 

What one finds , in the G = SU(5) case, is that each of the four irreps 1, 24, 75 and 200, into which (24 ® 24) s 
decomposes (A.l), contains exactly one standard model singlet $ w (the 4 eigenvalues of c ut are nondegenerate: 0, 5, 



,0 



8, 12). Therefore, in 577 (5), specifying the irrep of H determines its (standard model singlet) vev and subsequently 
its contributions to (A.8) up to an overall factor; in particular, the relative contributions between terms G^G^" in 
(A.8) are completely determined. 



For the 5O(10) unified gauge group, however, only the 1 and 54 irreps contain exactly one standard model singlet; 
the 210 contains three, and the 770 contains four linearly independent standard model singlets (the 9 eigenvalues of 
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group 


C* 2 (n) 


C* 2 (r 2 ) 


C 2 (r 3 ) 


SU(5) 


C* 2 (24) = 5 


C 2 (75) = 8 


C 2 (200) = 12 


50(10) 


C 2 (54) = 10 


C 2 (210) = 12 


C 2 (770) = 18 


SU(4) 


C* 2 (15) =4 


C 2 (84) = 20 




5(7(2) 


C* 2 (3) =2 


C 2 (5) =6 





TABLE V: Quadratic Casimir invariants C 2 (r) for some irreps r of several groups of interest, in the conventions described 
above. C 2 (l) = for the singlet 1 of any group. 



c u t are partly degenerate: 0, 10, 3 x 12, 4 x 18). These statements hold for both the normal (Georgi-Glashow) and 
the flipped embedding of G321 into 5O(10), and, in fact, it is found that the set of standard model singlets is the 
same in either case. 

In order to distinguish these linearly independent G321 -singlets within the 210 or the 770 of 5O(10), one can specify 
their transformation properties under subgroups of 5O(10) that contain the standard model G321 (the normal or the 
flipped G321, respectively). We consider here two such subgroups, namely the two maximal subgroups of 5O(10) 
consistent with the charge assignments for the standard model fermions; these are (a) the extended Georgi-Glashow 
subgroup SU(5) x U(l) x C 5O(10), and (b) the Pati- Salam subgroup SU(4) C x SU{2) L x SU(2) R . T o achi eve this 
practically, the method of the Casimir operator (A. 14 1 is employed again; however, the sum over c in (A. 14) is now 



restricted to the generators of the desired subgroup, or, more precisely, to the generators of each simple factor of the 
subgroup, in order to find linear combinations of the <& w that transform in irreps of the subgroup (factors). It turns 
out that either one of the subgroups (a) or (b) can resolve the degeneracy of the vectors <fr w transforming within the 
210 or the 770 of 50(10), see the second column of Tables VII and |VIIl] for cases (a) and (b), respectively. 

So, for each of the irreps r on the right-hand-sides of (A.l) in the 5(7(5) as well as the 50(10) case, we have 
found a basis {&( r )t} t i such that any standard model singlet tfi ab in the irrep r can be written as a linear combination 
4> ah = Ylt c t&(r)t with real coefficients Ct, the index t runs over t — 1 only, except for 50(10) in the cases of r = 210 
(where t = 1, 2, 3) or r = 770 (t — 1, 2, 3, 4), for which the 3?( r )t can be given definite transformation properties 
under (maximal) subgroups of 50(10), as described above. To be specific, we normalize the overall magnitude of 



each $ 



(r)t 



such that£«6*$ t *& t 



= 1 



resulting in ^)t^(r')t' = S r 



■'St., 



Also, it turns out that 3>?A t = whenever 



a and b belong to two subsets of generators with different transformation properties under G321, see (A. 9), ( |A.1Q 1 or 



(|A.11|), and ^f^\ t is proportional to the identity matrix when a and b are restricted to one set of generators. Therefore, 



nonzero off-diagonal elements ^(jJw 7^ can possibly occur only in the 50(10) case when a and b correspond to the two 



degenerate G32i-singlet generators III, V or III', V . In particular, $"*) t is diagonal in a, b on the set a,b — 1, . . . , 12 
and is proportional to the identity matrix on each simple factor of the standard model a, b = 1,...,8 (SU(3)c), 
a, b = 9, 10, 11 (SU(2)l) and a, b — 12 (U(1)y)', one could therefore arbitrarily, as an invariant convention, for each r 
and t choose the sign of 3?( r )t such that, e.g., the nonzero diagonal entry $?pw with smallest a = b is positive; however, 
to save some writing, we reverse this convention for the flipped embedding case in rows 5 and 8 of Table |VIII| 

We are now ready to write down the matrices obtained in this way, explicitly. Since, as just described, each 

matrix &( r )t contains only a few independent entries, they can be written down in an economical way (1^ denotes the 
identity matrix of size d x d): 



ab 



for 517(5) : $ Qb 



for 50(10) : $ 



ab 



V 



III 





iii.v Q 



ab 



(A.15) 



(A.16) 
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For the 5(7(5) case, the numerical entries of each matrix 3>( r )t = $m (A.15) are listed in Table VI (in this case, 



t = 1 only, see abov e) . Fo r both the normal and the flipped embedding of G321 C 50(10), Table 



entries of each $ 



(r)t 



VII lists the 



in' 



(|A.16| (for the flipped embedding, 3>( r ) t 



^(r)t an< ^ ®(r)t V snou ld be used from Table 



VII 



instead of the unprimed etc.), where the G32i-singlets (labeled by t) within each 5O(10)-irrep r are further 
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classified according to their transformation under the maximal subgroup 5(7(5) x U(l)x C 50(10), see above; Table 
is similar, but with the G32i-singlets &'r r \ t classified according to their transformation properties under the 
maximal subgroup SU{A)c x SU(2)l x SU(2)r C 50(10) (of course, all G32i-singlets are 5f/(2)_c -singlets, since 



SU{2)l C G321). In all cases, the matrices ( |A.16 ) with specified transformation properties under any subgroup are 
identical for the normal and flipped embedding except for the $ /// , <£> y and <J> //7 ' y entries (which correspond to 
generators t 111 , t v that have the same quantum numbers under G321), and these entries are related by 



III'V 




(A.17) 



with the 2 x 2-matrix V from (A. 12). 



r 


W 


W 


M 


M 


1 


l/v 7 ^ 


1/^24 


l/v 7 ^ 


l/v 7 ^ 


24 


2/V63 


-3/V63 


-1/V63 


-l/2\/63 


75 


l/v 7 ^ 


3/V72 


-5/V72 


-1/V72 


200 


1/V168 


2/v / 168 


lO/v 7 !^ 


-2/V168 



TABLE VI: The standard model singlets <f>( r ) in each of the irreps r ( A. 1 ) of 5(7(5) in the explicit version (A.15I with the 
conventions described above. 



50(10) 


5(7(5) x (7(l)x 






®(r)i 




*U 


(r)i 


®(r)i 




gyVIII 


®(r)t 


<) t 


(r)t 


r = 1 


(1,0), t=l 


1 

3^5 


1 

3^5 


1 

3^5 


1 

3^5 


1 

3^5 





1 

3^5 


1 

3-i/S 


1 

3^5 


1 

3^5 


1 

3^5 





r = 54 


(24,0), t=l 


1 

730 


/ 3 

V 40 


1 

2 v / 30 


1 

4^30 





1 

2v / 5 


1 

4^30 


1 

✓30 


/ 3 
V 40 


1 

2^30 





1 

2^5 


r = 210 


(1,0), t = l 


1 


1 


1 


1 


_ 2 





_ 1 


1 


1 


19 


2 


/ 2 


2\/15 


2\/15 


2Vl6 




\/l5 


2\/15 


2 V / 15 


2 V '15 


10\/15 


5\/15 


Y 125 




(24,0), t = 2 


2 


_ 1 


1 


1 





_ 1 


1 


2 


1 


7 


_ 4 


13 




3^15 


V15 


3^15 


6^15 


\/l0 


e-v/iS 


3^15 


vtS 


1S\/15 


5^15 


isvTo 




(75,0), i = 3 


1 

6^3 


1 

2^3 


_ 5 

6^3 


_ 1 

6^3 








1 

6^3 


_ 1 

6^3 


_ 1 

2V3 


_ 1 

30^ 


_ 4 
5^3 


72 
15 


r = 770 


(1,0), t=l 


1 


1 


1 


1 


/ 32 





1 


1 


1 


77 


/ 8 


/ 3 


6^10 


6^10 


6VT5 


6^10 


V 45 


3\/l0 


3^10 




30^10 


Y 1125 


Y 125 




(24,0), i = 2 


1 


1 


1 


1 





j~T 


/ 


/ 7 


/ 7 


101 


/ 128 


97 




3V210 


2 v / 2l0 




12V210 


V 20 


V 4320 


V 270 


V 120 


30V210 


Y 2625 


SOv^ 




(75,0), t = 3 


1 

6-V/6 


1 

2^6 


5 


1 

6\/6 








1 

3^6 


1 

3^6 


1 


1 

30^ 


Y 75 


1 

15 




(200,0), t = A 


1 

2^42 


1 

v/42 


5 

^42 


1 

s/42 

















1 


/ 96 

V 1 75 


2 

5^7 



TABLE VII: The standard model singlets $( r )j in each of the irreps r ( A.l I of 5O(10), classified according to their transformation 
properties under the 5(7(5) x U(l)x C 5O(10) maximal subgroup (second column; note, all standard model singlets $/ r ) t 



also happen to be (7 (l);y-singlets), in the explicit version (A. 16 1 with the conventions described above. The entries of the 
matrices *&( r )t (A. 16 1 agree between the normal and the flipped embedding of G321 C 5O(10), except for the &fr)ti ^(r)t an d 



&*1L' V entries, where the primed ones 



. 1 m j ,. .. „ ,., [„ ,„„ 1 ..„. v*(r)t i etc., see last three columns) should be used instead for the flipped embedding, 

cf. also | A.17| |. (For both the normal and the flipped embedding, the 5(7(5) here denotes the 5(7(5) C 5O(10) containing 
SU(3) C x SU(2) L .) 



Furthermore, for each 5O(10)-irrep r separately, the vevs $( r )t, classi fied a ccording to 5(7(5) x U(l)x (see Table 



VII I, and the vevs ^'( r )f> classified according to SU{4)c x SU(2)r (Table VIII), are linear combinations of each other, 



in tact, they are related by orthogonal transformations; when their overall signs are chosen as shown in the tables, 
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50(10) 


5(7(4) c x SU(2) R 


*(r)i 










(r)t 


®(r)t 


folVII 


&VII1 


«////' 

9 (.r)t 


<b' V ' 


(r)t 


r = 1 


(1,1), t = l 


1 

o V o 


1 

JV3 


1 

,5 V 


1 

ovo 


1 

375 





1 

ovo 


1 

OVO 


1 

ovo 


1 

ovo 


1 

375 





r = 54 


fl 1) £ = 1 


1 


/ 3 

V *o 


1 

2730 


1 

4^30 


o 


1 

275 


1 

4730 


1 

75o 


/ 3 

Y 40 


1 

2730 





1 

275 


r = 210 


(1,1), * = 1 
(15, 1) , t = 2 




1 

372 


i 

76 




/ 3 

_ 275 
15 






/ 2 

Y 75 
_ 72 

5 


1 

5 

_ 2 
573 









_ 1 
372 


1 

76 




/ 3 
Y 50 

_ 272 

15 


/ 2 

Y 75 

_72 
5 


1 

5 

_ 2 
573 




(15,3) , t = 3 








2 


1 


2 


1 

5vb 


1 








2 


2 


1 

5 v6 


r = 770 


(1,1), t = l 


1 


/_5_ 


19 


/ 5 


/ 32 


/ 3 


/ 5 


1 


/_5_ 


19 


/ 32 


~\/ ~ 


3\/lO 




30%/TO 










3710 




307l0 








(1,5), t = 2 








76 
5 





[Y 

Y 75 


2 
5 








1 

76 


76 

5 


Y 75 


2 
5 




(15,3) , t = 3 








272 
5 


1 

672 


272 
5 


1 

573 


1 

672 








272 
5 


272 
5 


1 

573 




(84, 1) , t = 4 


i 

6710 





/ 128 
V 1125 





/ 32 
Y 126 


8 

5715 





pi 

Y 45 





/ 128 
Y 1125 


/ 32 
Y 125 


8 

5715 



TABLE VIII: Same as Table VII but with the standard model singlets &t r ) t classified according to their transformation 
properties under the SU(4)c X SU (2)l X SU(2)r C 5O(10) maximal subgroup (second column; all entries are SU(2)l singlets, 
since SU(2) L C G321). (The 5(7(4) c here denotes the 5(7(4) factor in 5*7(4) x 5(7(2) x 5(7(2) C 5O(10) with 5(7(4) D 5(7(3) c , 
where one of the two 5(7(2) factors is SU(2)l.) 



then, for both the normal and the flipped embedding, the relation between the 3>( r )t and the &'t r u is: 



(1)1 



$ (1 



(1)1 



(54)1 



(54)1 



/*(210)1 
^(210)2 
\ < ^ ) (210)3/ 



^(770)1^ 




^(770)2 




V (770)3 




^(770 J 





V 



20 
3 

275 

3 
5 



3^ 
175 

4 
/35 



7105 

16 
572T 



2 
5 

8_ 
15 



3 
5 

_ 1 

2 

1 

273 
1 

715 







$(210)l\ 


1 

76 




$(210)2 


^73/ 




\*(210)3/ 







/ 35 




($(770)0 


1 

77 




$(770)2 


/? 
V 7 




$(770)3 


/IT 




V$(770) 4 / 



35 



(A.18) 



By classifying the standard model singlets $( r )t, in the SO(10) case, into irreps under (maximal) subgroups of 
5O(10) as in Tables VII and VIII we do not mean to imply a grand unified symmetry breaking scenario where 5O(10) 
is broken to the standard model G321 necessarily via some intermediate gauge group, although this classification is 
well suited for such a scenario, see the end of sectio n |VI[ Moreover, such a classification in particular can serve as a 



parametrization for standard model singlets (Hf °) (A 



transforming in irreps r,, in terms the of the basis vectors 



ab\ 



t 



Vi4>: 



ab 



with real scalars Vu\ t , and where, e.g. for = 210, 



i a b 



$$ 10)1 cos 9 l + ^ 10)2 sin d x cos 9 2 + $ ab 



(210)3 



sm U\ sm 02 , 



(A.19) 



(A.20) 



such that 4>f b (f>f b = 1 and J^t v fi)t ~ v 1 f° r a ^ then the kinetic terms of the Vi and v^ t field degrees of freedom, before 



assuming vevs, are canonical (A.3) C = \ (d^H. 



Tab 



ab 



h {dfiVi) 2 + . . .. (The two parametrizations 



in (A.19) differ only if Hi is a 211) or 770 of 5O(10), since in all other cases the sum merely runs over t = 1.) 



We assume, for definiteness and simplicity (see Sections [TT| and VI), one-step breaking of the grand unified gauge 
group G to the standard model G321 at the unification scale Mx- Below Mx, all Higgs multiplets Hi responsible for 
grand unified symmetry breaking assume nonzero vevs (A.19), and in particular give masses to the non— G32 1 g^uige 



bosons (henceforth called "superheavy" gauge bosons, although not all of them can get mass from the Higgses ( A.l ), 



24 



see below), see (A.3) 



C = l -g 2 G AlA b » E Tr (- [t a G , (H,)} [t b Gl (HA]) ee \ml h A«A b » 



(A.21) 



each Higgs Hi contributes independently to the gauge boson squared mass matrix m 2 ab . For an SU(5) grand unified 
group, each of the 12 superheavy gauge bosons acquires equal mass: 



51/(5) 



V- C 2 {r l ) 2 2 

2^ 9 G Vi Sab 



Mi 5 ab for a,b= 13, 



,24 



(A.22) 



where C^r,) is the quadratic Casimir invariant of irrep Ti (cf. Table [V] for the irreps ( |A.1[ )), and Vi the vev of Hi 
(cf. ( |A.19 ) and below). In SO(10), the situation is more complicated: Higgs multiplets in any of the irreps ( A.l I of 
interest for the dimension-5 operators (A.7) fail to give mass to all of the superheavy gauge boson s. In p arti cular, the 



non-G32i gauge boson G v ^ v (or GY V for the flipped emb eddin g) belonging to the set V (or V; see (A. 10) or (A. 11)) of 
generators is always left massless by Higgses in irreps (A.l I, provided only that their vevs are G32i-smglets, ct. also 



the second column of Table [VII| Also, some of the other superheavy gauge bosons might remain massless, and, at any 
rate, do not receive equal masses. Higgses Hi in representations other than those on the right hand side of ( A.l ) are 



needed to give mass to all of the 33 non-G*32i gauge bosons in the 5 *0(10) case, and they cannot occur in dimension-5 
operators (A.7). Independent of the individual directions 4>f b (A. 19) of the vevs (Hi) = Vi4>i, the averaged superheavy 



gauge boson squared mass is: 



SO(10): m\ b - 



EC 2. ( f*i ) 2 2 2 
33 9g v i s ab = M gb 



5 a b for a,b = 13, 



,45 



(A.23) 



where the sum runs over all Higgs multiplets at the grand unification scale, whether or not they occur in (A.7). 



Grand unification asserts that the masses of the superheavy gauge bosons are related to the scale Mx of gauge 
coupling unification; namely, their masses have to be somewhere around Mx, such that the renormalization group 
equations of the standard model apply below the scale Mx , at which the gauge couplings unify and above which the 
unified gauge coupling evolves according to the /3-function of the unified gauge theory. For definitcness in obtaining 
numerical values, we assert for the analysis in the main text, that the superheavy gauge boson masses (or the averaged 
superheavy gauge boson mass (A.23) in the 50(10) case) have to exactly equal the unification scale: 



M, 



gb 



M 



x 



(A.24) 



In the SU(5) case, any of the Higgses in (A.l) can give mass to all of the superheavy gauge bosons; so we assert, 
again for definitcness (see Sections [Tl| and III[) , that no Higgs multiplets in irreps other than (A.l) are present in the 
theory (nor, for that matter, in the sum (A.22)). In the SO(10) analysis, however, we assert that the Higgs multiplets 
in irreps (A.l) account for some fraction 1//, e.g. half, of the average gauge boson squared mass: 



M 



x 



M 



\ 



d{G) - 12 



{i | r-j in dXT}} 



9GV 



\ 



f 



d{G) - 12 



(A.25) 
(A.26) 



{i | Ti in fATTl} 



with / = 1 (resp. / = 2) for SU(b) (resp. 5O(10)), and d(G) — 12 = 12 (rcsp. = 33) the number of superheavy gauge 
bosons. The form (A.26) arises when the vevs = XiV are to obey some given ratio x± : x% : ... , with an open overall 



scale v, which will be useful later. 



Now, plugging the Higgs vevs (A. 19 1 back into the Lagrangian (A.8) and adding the gauge boson kinetic term 
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( A. 3 1, one obtains: 
C 



1 



G*G aiJ - v 




y-v su(3) 



afif 
SU(2) 



r ^ r U(l) 



-(l + e 3 )F« v F: 



SU(3) 



(1 + e 2 )K,F. 



SU{2) 



-(l + e 1 )F llv F^ 1 



(A.27) 



Here, we have used the forms (A. 15 1, (A. 16 1 of the standard model singlets &u)t on the subspace I, II, III (or 777') 
of G32i-generators and have omitted the superheavy gauge bosons, assuming the Lagrangian (A.27) to be valid below 
the unification scale Mx where they have acquired mass and are integrated out; the last line defines what we mean 
by e s , s = 1,2,3, namely the corrections to the gauge boson kinetic terms originating from the Higgs multiplets in 
the dimension- 5 operators assuming nonzero vevs {A 



e * = E E ( for s = !. 2 ' 3 )» where 5 i 



cot _ 



in 

(i)* 



:(i)t _ ^11 
°2 = ~^(i)t 



(A.28) 



the can be read off from Tables 



VI 



VII 



and 



and 



VIII 



Table |ITf in t 



for the S77(5) and the 50(10) cases (use 

le main text gives them for SU (5) explicitly (in this case, t = 1 



for the 



flipped embedding G 32 i C 50(10)), 
only). 

As shown in Section III after rescaling — > (1 + e^^F^ and g s — > (1 + e s ) _1 / 2 .g s to the observed low-energy 
gauge field strengths ana gauge couplings by amounts which are different for each of the standard model gauge group 
factors s = 1, 2, 3, the condition for gauge coupling unification at the scale M x into the unified group reads, in terms 
of the observed (running) gauge couplings a s = a s (fj,) = g s (fi) 2 /4tt of the theory below Mx- 



(1 + ei)ai(Mjf) = (1 + e 2 )a 2 (A/x) = (1 + f- 3 )a 3 (M x ) 



9h 

47T 



(A.29) 



where go = gc{Mx) is the gauge coupling of the unified theory at the scale fi = Mx- The functions a s (fi) are fixed 
through their low-energy measurements (e.g., at the scale fi — mz of the Z mass) and their renormalization group 
evolution; in the non-supersymmetric standard model to one-loop order with one standard model Higgs doublet: 



a s {fj) a s {m z ) 



2tt mz 



with /^-function coefficients 



h =10' 62 



19 

"6 



and initial values [50] 

ai (m z ) = 0.016887 ±0.00040, a 2 {m z ) = 0.03322 ± 0.00025 , a 3 {m z ) = 0.1176 ± 0.005 



(A.30) 



(A.31) 



(A.32) 



for the couplings at the scale mz = 91.1876 GeV. For two- loop evolution, which we use in Section III A for comparison 
to one- loop, see [17] . 

Under the perspective of this Appendix, unification is the numerical requirement that the unification scale M x , the 



unified gauge coupling go, the Wilson coefficients a and the Higgs vevs vf — J2t v (i)t sa tisfy equations (A.29) 



and (A.25) simultaneously, with the given running coupling functions a s (n) of the standard model, e.g. (A.30 



to one loop. In the main text, we further demand that the values of these quantities are "natural", so that p 
sensible unification can be claimed, and we exhibit such models. 



A.28 

{OS 



rysically 



2G 



We outline here a procedure to solve these equations simultaneously, after having chosen the unified gauge group 
(SU(S) or 5*0(10)) and the Higgs content Hi (the multiplets that may occur in the dimension-5 operators (A. 7)): 
Fix any ratio X\ : x% : . . . between the Higgs vevs Vi — xiv (without loss of generality, J^. xf — 1 with, e.g., spherical 
coordinates); for Higgs multiplets Hi in a 210 or a 770 of S'O(IO), also fix th e direction of the vev (Hi) by fixing the 
ratio U(i)i ■ U(i)2 ■ ■ ■ ■ between the v^t — U(i)tVi (with y%\ t — 1), see (A. 19). Also fix the ratio Z\ : z% : . . . between 
the Wilson coefficients a = ZiC (with J^i z f = !)■ With these inputs, the ratio between the e s (s = 1, 2, 3) in (A.28) 
is completely fixed: 



cv 

Mpi 



iV{i)t 



(A.33) 



A key observation is now that any such given ratio t\ : £2 : £3, along with the given functions a s (n), determines Mj 
and the values e s uniquely by solving the two left equalities in (A.29); this can be seen analytically to one-loop order 
when (A.30) is plugged back into (A.29) (the linear system of equations for InMx has nonvanishing determinant 
for almost any given ratio £1 : £2 ■ £3), but also holds at two-loop order. Then, the unified coupling ac — 9g/47t is 
determined by (A.29) as well, and so is the overall scale v of the vevs via (A. 26 ). Finally, with a choice ^ of the 
Planck scale Mpi, the required Wilson coefficients q = z^c can be computed from ( A.33| ). Note, that this last step is 
the only place where the choice of the Planck scale Mpi = 1.2 x 10 19 GeV/£ comes in: the smaller a Planck scale one 
chooses, the smaller the Wilson coefficients Ci have to be in order to achieve unification (at the same g/mt, Zi); in 
fact, the necessary Wilson coefficients are inversely proportional to the choice of £ in ([2]). 
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